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Prólogo de Editor 

Esta obra surge del material de clase elaborado durante los años 2016 y 2017 
para el dictado de las materias: Aplicaciones Industriales del Ultrasonido y 
Ultrasonido para Uso Médico, de las carreras de Ingeniería Electrónica y 
Bioingeniería respectivamente, dictadas en el Departamento de Electrónica y 
Automática, Facultad de Ingeniería, Universidad Nacional de San Juan. 


Se trata de transparencias que resumen los contenidos teóricos y prácticos del 
ultrasonido. Estos contenidos han sido elaborados a partir de la obra 
"Ultrasonidos para Ingenieros y Estudiantes de Ingeniería", bibliografía general de 
ambas cátedras, que se pueden obtenerse en forma gratuita en el sitio web: 
http://dea.unsj.edu.ar/ultrasonido (Biblioteca de cátedra). 


Si bien el material ha sido elaborado de forma tal que resulte claro, el lector 
deberá apoyarse en la obra arriba citada para una mejor y rápida comprensión. 


La obra está compuesta de 7 temas que se resumen a continuación: 


Tema l: se hace una revisión histórica de los avances y aplicaciones del 
ultrasonido, mostrando en forma muy sintética los temas que se desarrollarán en 
los capítulos posteriores. Se tratan cuestiones básicas acerca de los tipos de 
ondas, de la reflexión y trasmisión del ultrasonido en diferentes medios materiales 
y biológicos. Surge el concepto de impedancia acústica y se repasan las bases 
físicas del efecto Doppler. 


Tema Il: se analizan los diferentes métodos de registro usados, tanto en las áreas 
de los ensayos no destructivos como en medicina. Se muestran los esquemas 
básicos de sistemas de emisión y recepción para las dos áreas, sentándose las 
bases de la ecografía industrial y de la ecografía médica. En el área de los 
ensayos no destructivos se analizan con detalle, dos sistemas de probada 
eficacia. 


Tema lll: se encuentra dividido en dos partes; la primera que se refiere a los 
transductores de ultrasonido de bajo nivel de potencia, se analizan todas sus 
características. Se hace una correspondencia física entre las lentes ópticas y las 
acústicas. En la segunda parte, se profundiza el concepto de array ultrasónico, se 
realiza una revisión de los métodos para el control de la focalización, tanto en 
emisión como en recepción. Se muestra cómo a partir de la simulación, se puede 
estimar y modificar las características del campo acústico en las cercanías del 
array. 


Tema IV: se aborda el modelado de la señal de eco. Se presentan un número 
importante de variantes, tanto en el dominio del tiempo como en el de la 
frecuencia y se analizan las bases de la reducción de datos. Se da suma 
importancia al ambiente de simulación. 


Tema V: trata el análisis de los datos dados por los fabricantes de transductores 
para aire, palpadores para ensayos no destructivos y para medicina. A lo largo de 
una serie de ejemplos, el lector aprende a estimar las características mínimas que 
debe cumplir el sistema de ultrasonido, para una aplicación específica. 


Tema VI: se sientan las bases del efecto Doppler aplicado a la medicinal. A lo 
largo del mismo, se analizan los métodos y las arquitecturas para la generación y 
recepción de Doppler continuo y pulsado y, cómo se superpone la información de 
velocidad con la ecografía de referencia. 


Tema VII: se presenta a los transductores ultrasónicos de potencia. Se parte del 
estudio de las características de los materiales cerámicos y metálicos que 
componen un transductor de este tipo. En función de las analogías mecano- 
acústicas-eléctricas se analizan los modelos clásicos y se sientan las bases para 
el diseño. Al final del capítulo se aborda el tema del sonar. 


En toda la obra se hace especial hincapié en el empleo de la simulación como 
herramienta idónea, tanto en la faz de verificación de la teoría como en la del 
diseño. El software que se emplea debe tomarse como referencia, ya que no es 
único. 

El lector interesado en esta área del conocimiento puede visitar la página web 


anteriormente citada donde encontrará gran cantidad y variedad de material 
bibliográfico, o escribirme a: avecaGOinaut.unsj.edu.ar 


Angel C. Veca 
San Juan, Argentina, marzo de 2018 
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Que es el Ultrasonido (1) 


El ultrasonido son ondas acústicas, es decir, perturbaciones de tipo elástico, que se propagan 
por todos los medios materiales mediante el movimiento armónico de sus moléculas. Se 
llama así, porque su frecuencia de oscilación está por encima de la frecuencia de las ondas 
audibles, esto es, frecuencias superiores a 20kHz. El límite superior no puede determinarse 
pues depende de la tecnología. En este contexto, actualmente se está trabajando a 
frecuencias superiores a los 100 MHz, (microscopía ultrasónica). 


pop 


Frecuencia 0,0...1Hza14,99Hz 15Hza1999KHz 20 KHz en adelante 


Que es el Ultrasonido (2) 


Ciertas aplicaciones utilizan la potencia de las ondas ultrasónicas para transformar el medio 
al que se aplican. 


Rango >0,1 W/cm? 0,1 W/cm? a KW/cm? 


Tal como ocurre en medicina donde son usados con fines terapéuticos como lo son la 
hipertermia, litotricia, etc 


. En la industria, donde su uso más frecuente es para activar procesos de naturaleza física 
como filtros, aglomerantes, destrucción de burbujas en líneas industriales de 
embotellamiento de bebidas gaseosas. En la industria química para el control de la 
aceleración de reacciones. Así como también limpieza y soldadura de piezas, entre otros 


Que es el Ultrasonido (4) 


Otro campo de aplicación consiste en interpretar los cambios producidos en las ondas 
ultrasónicas a su paso por el material, permitiendo deducir las características mecánicas y/o 
geométricas de éste, para determinar las propiedades físicas de la producción y determinar la 
existencia de defectos internos, tales como poros, grietas, calidad de las uniones, etc. Este es el 
caso de los ENSAYOS NO DESTRUCTIVOS (END). 
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Es un campo ilimitado de aplicaciones ya que, actualmente se utiliza para el control de 
calidad en las industrias aeroespacial, nuclear, siderúrgica, naval, automotriz y de 
construcciones civiles. 


Que es el Ultrasonido (5) 


Como se desprende de lo anterior y de lo que se verá en este curso, el ultrasonido es una 
herramienta muy simple de emplear y sobre todo muy poderosa, ya que, permitirá "ver" 
donde le es imposible a la visión (normal o artificial) y no producir cambios en los 
organismos vivientes como sucede con los rayos X, dado que se trata de una radiación no 


i¡onizante. 


Revisión Histórica (1) 


1793 - L. Spallanzani estudió el sorteo de obstáculos de los murciélagos en la oscuridad. 
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Revisión Histórica (2) 


- F. Galton. A él se le debió entre otras, el desarrollo del 1er em | y 
las primeras pruebas con seres humanos y con animales para ys de 
, afrecuencias por debajo y por arriba de 


Se sugiere el empleo del ult la de ¡ s después del 
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Revisión Histórica (3) 


1929 - Se emplea el ultrasonido para detectar presencia de heterogeneidades en materiales, 
dando el primer paso al método de transparencia y a los ensayos no destructivos con 
ultrasonidos (END US) 


1939 - la tecnología electrónica da un paso fundamental, nace el Radar, aportando 
innumerables mejoras al Sonar. 

1947 - Se sientan las bases de la ecografía médica. D. H. Howy construye partes de lo que 
seria el primer ecógrafo bidimensional para uso médico. El órgano a estudiar debía 
sumergirse en un tanque con agua. 


Revisión Histórica (4) 


1948 - Se aplica el principio del sonar para la detección de defectos en materiales, 
iniciando la técnica pulso-eco en ensayos no destructivos. Concepto de A-scan 
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Revisión Histórica (5) 


1955 -— Se descubren las propiedades piezoeléctricas del PZT. 


1957 — Se construyó un escáner de contacto bidimensional, evitando así la técnica de 
inmersión. 


1959 -— Se sienta las bases de la velocimetría Doppler en ecografía 


Registro de flujo 
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Revisión Histórica (6) 


1971 - Se desarrolla el primer array lineal. 
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1983 - Se introduce al mercado el primer Equipo de Doppler a Color, que permitió 
visualizar en tiempo real y en color el flujo sanguíneo 
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Ultrasonido en la Naturaleza 


Espectro ultrasónico en la naturaleza 
Rango: 18 kHz a 100 kHz 
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Síntesis de las Bases de la Materia (1) 


Concepto de Pulso-eco (1) 


o e Fuente de emisión direccional 
Fuente de emisión omnidireccional 


e 


Proceso de emisión y recepción: Pulso-eco 


d | d 
V=>—  Sicesla velocidad del ¿_“% 
sonido en el medio Í 
Nos interesa conocer la d 
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C 
Tiempo de transito: tiempo de ida + de vuelta tf = 2d 
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Síntesis de las Bases de la Materia (2) 


Concepto de Pulso-eco (2) 


Modelo elemental del proceso de emisión-recepción 


Secuencia emisión-recepción 


Pulso de 
emisión 


| 
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Síntesis de las Bases de la Materia (3) 


Propiedades del Medio (1) 
Concepto de Impedancia Acústica (1) 


os lie fe 
AN o > e a es 
pa jo o Lo [Je E 
Impedancia eléctrica Impedancia acústica 
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Unidades — |(Z,)= ] 
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Síntesis de las Bases de la Materia (4) 
Propiedades del Medio (1) 


Concepto de Impedancia Acústica (2) 


Medio | plKgm*]| ¿Em m] | ZlKgm*s"] 
+10? Cr m7 Es a +1 p0 
+10? 107 


a CO IC O A 
Agua 1.00 L3 13 

Resinas y Plásicos | EL 
Aluminio 2.71 6.32 2,2-3,2 17.1 


Bronce 7.4-8.9 | 4.4-4.98 32.5 — 44.5 


Glicerina 126 | 192 [192 | 24 
Perspex (Acrílco) 1.18 2:13 1.43 32 
Tolueno 0.87 | 133 [MN 09-12 


Aceites 0.9 1.2-1.7 1.2-1.7 
Maderas 0.6-1 1.4 0.5-3 
Vidrios 2.2-3.9 2.56 — 3.42 2.56-3.52 12.3= 15,1 


Tejidos humanos 1.06 1.54 1.63 


_pP 
La = 7 Reemplazando Pyv | Za=p.cC 
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Síntesis de las Bases de la Materia (5) 


Propiedades del Medio (2) 
Coeficientes de Reflexión y Transmisión (1) 


medio 1 []  medio2 ¿Qué porcentaje de la presión emitida se refleja retornando 
A a la fuente y que porcentaje atraviesa la interfaz? 
| | da > Dz 
e. R, Coeficiente de Reflexión — R=“=“2 1 
| ii Z,+Z, 
Onda transmitida 
Onda incidente a 
H ' . . 7 Zo 
T, Coeficiente de Transmisión “Y” = ———“— 
LL 


Z,, Impedancia Acústica del medio donde se propaga la 
onda incidente 

Z,, impedancia Acústica del medio donde se transmite la 
onda incidente 


Interfaz: Cambio de medio 


£L, 
Ejemplo de tres medios 
Dos interfaces > 


Síntesis de las Bases de la Materia (6) 
Propiedades del Medio (2) 


Coeficientes de Reflexión y Transmisión (2) 


Ejemplo 

Detección de bancos de peces. La impedancia acústica del tejido del pez es similar a la del 
agua. La pregunta es inmediata: ¿Cómo es posible la detección de peces? El sonar no detecta 
al pez en sí, lo que detecta es el aire en el interior de su vejiga natatoria. Para verificar esta 
aseveración, calculemos el coeficiente de reflexión de presión para una interfaz agua — aire, 
sabiendo que las respectivas impedancias acústicas son: 


_ o Kg 
dE 428-1,5-10* 
por lo que: » SAT E 0,99 
hs 428+1,5-10 
Zhieto =2,95:10* 
el 99% de la presión incidente se refleja en la interfaz 
Ejemplo 02 02 
Interfaz Agua-Músculo-Agua R, A 0,0625 R, pe e -0,0625 
la 2 0 
Z =1,7-10" 
musculo m?s 3,4 3 
T, =—=106 TL => =0,938 


Síntesis de las Bases de la Materia (7) 


Propiedades del Medio (3) 
Atenuación de la Energía Acústica (1) 


A medida que transcurre el tiempo, la energía acústica emitida se 
distribuye en una superficie mayor 


1 2 
; gal an 0d f 
a a P=P —€ 
5, | S, S S, d 


Po es la presión de emisión, d la distancia 
recorrida total, o. la constante de atenuación 
del medio y f la frecuencia de trabajo 


fl< A < 13 
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Síntesis de las Bases de la Materia (8) 
Propiedades del Medio (3) 


Atenuación de la Energía Acústica (2) 


Ejemplo de la atenuación en tejidos 
biológicos MELO o.[dB/cm] 


Valores típicos de ecografías: 
Abdominal: 2.25 a 3 MHz 


Ecocardiografía: 3.5 MHZ Pequeñas partes: 7.5 a 10 MHz 


Ecocardiografía pediátrica: 5MHZz al 


Síntesis de las Bases de la Materia (9) 
Propiedades del Medio (4) 


¿Que es lo que el murciélago no puede hacer? 


1-Detectar una mariposa detrás de una hoja 


Hoja 


hb - Ml ++ -) -) Mariposa vb E Ñ Y a —= Mariposa 


2-Detectar una mariposa a una distancia mayor de 5 metros 


al = Asta + A emnelta + Anottexión 
Asotal > 60dB 


3- Focalizar el haz de ultrasonido 


Tiene que mover la cabeza para orientar el haz de ultrasonido 7 


Síntesis de las Bases de la Materia (10) 
Naturaleza de las Ondas Ultrasónicas (1) 
Tipos de Ondas Ultrasónicas 


Ondas Longitudinales — Transversales - de Superficie 


As. 
=> - “XEAR> 


Simétricas — Anti simétricas 
Ondas de Lamb - Transversales - de Superficie 


Velocidad de propagación de las ondas Longitudinales 
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Síntesis de las Bases de la Materia (11) 


Naturaleza de las Ondas Ultrasónicas (1) 
Tipos de Ondas Ultrasónicas 


Ondas de Presión u Ondas P u Ondas Longitudinales, (L) 
Se transmiten en medios sólidos, gases y fluidos 


Ondas de Corte u Ondas S u Ondas Transversales, (T) 
Se transmiten sólo en medios Sólidos. 
Su velocidad de propagación es menor que el de las Ondas P 


Ondas de Superficie 
Se transmiten sólo por la superficie de sólidos 


Ondas de Lamb 


Cuando el espesor del material es inferior o del mismo orden de la longitud de onda de 
las ondas ultrasónicas 
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Síntesis de las Bases de la Materia (12) 


Naturaleza de las Ondas Ultrasónicas (2) 
Reflexión de Ondas Ultrasónicas (1) 


E | 


H: conocido 
H1, H2, L,L1 y L2 muy grandes 
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Síntesis de las Bases de la Materia (12) 


Naturaleza de las Ondas Ultrasónicas (2) 
Reflexión de Ondas Ultrasónicas (2) 


Oli | 


H: conocido 
H1, H2, L,L1 y L2 muy grandes 
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Sintesis de las Bases de la Materia (12) 
Naturaleza de las Ondas Ultrasónicas (2) 


Reflexión de Ondas Ultrasónicas (3) 


H1 Palpador de US 


H2 ] : 
¿? (04 


ÍÁí—_—_—__— A ----—_—_—— 
kK—áKáÁá- 11 + 
WA 


¿Cuál es el ángulo OÍ que tiene que formar el haz de ultrasonido para que se detecte 
un defecto esférico? 

Según la forma y posición del objeto, puede producirse difracción y/o cambio de 
modo, debiendo estudiar la técnica que sea más idónea. 


T 
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Síntesis de las Bases de la Materia (15) 


Reflexión y Transmisión de las Ondas Ultrasónicas (1) 


Ley de Snell Interesa conocer 0%, 


medio 1 medio 2 C 
j t A senQ, == sena, 
a o, senó, C, C, 


Rápido 
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Síntesis de las Bases de la Materia (16) 


Reflexión y Transmisión de las Ondas Ultrasónicas (2) 


Agua Aluminio Aluminio Agua 
| medio 2 | 


medio | . medio 2 


E 


Solido Solido 
medio 1 ; medio 2 
e a e 


Síntesis de las Bases de la Materia (17) 


Reflexión y Transmisión de las Ondas Ultrasónicas (3) 


Casos particulares 


Solo ondas L Solo ondas T Ondas L y T 


sen, , A send, , send, ; A send,., send, ; A send, ; 


Co Cr 


senO, , sen0,, send, ,  sen0,, 


Ley de Snell Generalizada a 


Síntesis de las Bases de la Materia (18) 


Concepto de Efecto Doppler (1) 


Concepto físico Casos: 
a: cial 1- Emisor y receptor moviéndose en la misma dirección y sentido 
- Mayor Menor 
q. a” 


SS o a 


— 2 — 
S: Emisor, R: Receptor, Vs: velocidad del emisor, Vm: velocidad del 
A medio, V¿: velocidad del receptor, fo: frecuencia de emisión 


2- Emisor y receptor moviéndose en la misma dirección y sentido 


opuesto 

5 Vs Mb O R Ú 
03 J))))) 1.548) 
sx Ss Y sx 


- 
a e a 


31 


Síntesis de las Bases de la Materia (19) 


Concepto de Efecto Doppler (2) 


3- Emisor en reposo y receptor moviéndose en la misma dirección y ambos sentidos 


A S a . e R e 
MX Ve > 
1 , F R R 4 
0: ))))))).0—, 
ES 
tf V 
E i 
2 fecos O 
ES 


2 


Síntesis de de la Materia (20) 
Concepto de Efecto Doppler (3) 


. . ] ¡ ; ¡ 
Universo del objeto Universo del sistema 
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Aplicaciones Industriales del Ultrasonido 
Ultrasonido para Uso Médico 


Prof. Dr. Angel C. Veca 
Marzo 2018 


Tema ll 
Técnicas de Emisión y Recepción 


TEMA 


Técnicas de Ultrasonido continuo 

Técnicas de Ultrasonido Pulsado 

Zona ciega 

Método del Umbral 

Registros típicos 

Pantallas y Modos 

A-scan. Casos Típicos 

Sistemas para A-scan 

Diagrama en Bloques de un Sistema para Ecográfica Médica 
Relación entre A-scan, B-scan y Modo M 

Control Atenuación Distancia 

Referencias 35 


Técnicas de Ultrasonido Continuas 


Esquema implementación Método de 
Resonancia 


Resonancia 


Osciloscopio 


Oscilador senoidal de 


1 


; 
Ep frecuencia variable 


Oscilador senoidal de 
frecuencia variable 


Ensayo por transparencia y 
E ES acoplamiento con agua 


Lt. 
A 
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Técnicas de Ultrasonido Pulsado 


Robótica END-MEDICINA 


Ñ hm) 


Transmisión o transparencia (Sólo en END) 


Son 


Propiedades del Medio 


Coeficiente de Reflexión 


Impedancia Acústica Z= pC 


R qn L ¿0% L; 
La Es 
Coeficiente de Transmisión 
__2Z, 
Li + Lo 
¡| -mamp> Firma del Medio 


¡A A a A A O A a a A AA A 


| Eco de Fondo (EF) 
! a 
--4—-4— Ap A-Scan - Traza 


Amplitud - distancia 


xr (0) =xp(0) * h(t) = Rypx(0)Ó(t— Tp) + Apx()Ó(t — To Tr) +... 38 


Técnica de Pulso-Eco 


Esquema Básico Técnica de Pulso-Eco 


Emisor de 
Ultrasonidos 


Real 


Ideal 


topped 


reference 


| tert EZ rignt 


843.000 us ————1,09300 as 1.34300 as 
50.0 us/div 


nstante previo a la Comienzo de la 
interferencia interferencia 


2.29000 ns s 4.79000 ms 
500 us/div 


TINEBASE 


— 


l- delay 
LL 1.09400 ms] 


| reference J reference 


Ii tert MER right | ll | l I ! I | | nn ATA vet GEES right 
DN INN | 


fp stopped 


| Mindow NN AU | Il pl | | 
11 O on A | | PU a —— 


| ól ==) ad 


l O ] LEA ¡E 1 EA 
94.0000 us 1,09400 ms 2.09400 ns 94.0000 1.09400 as 2.09400 as 
200 us/div 200 us/div 


Evolución de la interferencia 40 


Proceso de Formación del Eco y Zona Ciega en END 


Inspecciones Automatizadas con Incidencia Normal e Inmersión 


L, >> £; El vaz 


Eco de Interfaz (El) 


Eco del Defecto (ED) 
DZ 


MATA 


Eco de Fondo (EF) 


A-Scan resultante 


41 


Pulso-Eco y Detección por Umbral 


Método Pulso-Eco por umbral Arquitectura del método 


Lógica 


de 
Control 


Traza amplificada Salida del detector amplificada 
Salida del detector 


Ll Í 
1.60000 as 
500 us/div 


a b 


Pulso de disparo y envolvente de la traza 
formada por dos ecos 
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Zona ciega saturada y eco 


reference 


left MA rignt 


Hindou 


2.20000 ns ————— 4,70000 as 
500 us/div 


44 


he stopped TIMEBASE 


Pulso de disparo 


Inhibición zona 
muerta 


4.14000 as  9,14000 as 
1.00 ms/div 


Pulso de disparo, inhibición del receptor y 
traza compuesta por un eco 
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Método de Pulso Transparencia 


Arquitectura del Método 


Ejemplo del Registro obtenido 


Indicación 


Mp cr coronel E 46 
A A A A A 


A E: E 


Sistemas de Representación 


Modo A o Pantalla tipo A o A-scan 
Amplitud de la Traza 


Modo B o Pantalla tipo B o B-scan 
Brillo o Intensidad de la Traza. Modulación de la amplitud 
Tiempo de transito y movimiento X-Y 


Modo C o Pantalla tipo C o C-scan 
Amplitud o Movimiento X-Y 


Modo D o D-scan 
Profundidad a la que se encuentra el máximo de la Traza en una Ventana 
Temporal pre-programada 


Modo M o o Pantalla tipo M 
Movimiento 
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Interpretación de Registros Típicos (1) 


Caso |: Medio sin impurezas 
ni defectos 150 


A 
E 
E A 


L/ 


0 500 1000 1500 2000 2500 3000 
Caso Il: Medio conteniendo 
una impureza e Tornillo 3cm del fondo a Tornillo 1.5cm del fondo 
50 
0 


Interpretación de Registros Típicos (2) 


Caso Ill: Medio conteniendo 
una impureza y una fisura 


Caso IV: Medio con una impureza o 
fisura inclinada respecto al haz 


o | 


Caso V: Medio con falta de 
paralelismo de las caras 


7 


oe 


Caso VI: Medio que contiene 
microporosidades 
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Interpretación de Registros Típicos (3) 


Ejemplos de Ecos en un medio altamente dispersivo 


500 1000 1500 Muestras SU 


Interpretación de Registros Típicos (4) 


Caso VII: Medio con superficies paralelas de reducido espesor 
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Esquema de un Sistema Analógico para A-scan 


SELECCION 
DE 
FILTROS 


E DEMO 
LIMITADOR | > T— DULADOR 


| 


PORTADORA 
MODULADA 
DEL ECO 


MCDULANTE 
DEL ECO 


GANANCIA 


BASE DE TIEMPO 


Y PANTALLA 
SINCRONISMO 


AH | Xx Yi 


2 


Ejemplos de Registros con Sistemas para A-scan 


Registro de una traza correspondiente al eco devuelto por un envase de 
puré de tomate 


PORTADORA 
MODULADA 


MODULANTE 
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Sistema SENDAS (1) 


Sistema para Ensayos No Destructivos De Arquitectura Segmentada 


Desarrollado en la década de los años 90' en el Instituto de Automática 
Industrial (IAl), Consejo Superior de Investigaciones Científicas (CSIC) — 
Madrid, España 


El 
Da ea 


DE ENTRA DA aro 


CAPTURA IMMDIATA DE DATOS 


CIAO ME 
ess en 


0dB 2120 de 
Daz Mhz 


FILTROS 
RF 


ENWOLWENTE 
FILTROS 
PA SABAJOS 


AL PALPADOR 


PULSER 
100 2600 Y 
05a2Mhz 

1316 PUL50O5 

FLA NCO + 


16 CURVAS 
DE 


512 PUNTOS 
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Sistema SENDAS (2) 


Módulos del Sistema SENDAS 


MODLIN 

Módulo de procesamiento analógico. Soportaba dos amplificadores logarítmico y/o 
lineal con un ancho de banda de 35 MHz. 

En ambos, la ganancia era programable (por soft) entre O dB y 120 dB. La salida de 
cualquiera de ellos podía ingresar al detector de envolvente para obtener la 
modulante. Ésta podía ser filtrada por medio de tres filtros pasa bajos (también 
programables por software). La elección de ésta dependía de la aplicación. 


MOCAD 

Módulo que implementaba la función CAD que podían generarse y grabarse hasta 
16 curvas diferentes y recuperarlas en cualquier momento. Cada curva tenía una 
longitud de 512 puntos. 


MODPUL 

Consistía en un pulser con tensión de excitación, frecuencia de repetición, número 
de pulsos y flanco de arranque programables: 100 a 600 Volt, 0.5 a 25 MHz, 1 a 16 
pulsos y flanco positivo o negativo, respectivamente. Tenía capacidad de manejar 
hasta 10 palpadores en forma cuasi simultánea. 


33 


Sistema SENDAS (3) 


MODAD 

Este módulo estaba compuesto por un conversor A/D con capacidad de muestreo de 80 
Mega y 8 bits de resolución. Debido a la gran velocidad de muestreo llevaba adicionado 
una memoria tipo FIFO para adecuar la transferencia de datos al procesador local. 


MODDI 

Este módulo se encargaba de realizar diferentes tipos de pre procesamiento de 
señales, por ello estaba compuesto por los siguientes sub módulos: 

EMI. Filtrado no lineal. 

DCV. Deconvolución. 

ENV. Extractor digital de envolvente. DKP. Detector monopico y multipico. 


La velocidad de procesamiento y transferencia de información entre los diferentes 
módulos era de 10 MS/Ss y todas las acciones de selección de parámetros, las realizaba 
el operador a través de la interfaz a usuario (Ver siguiente transparencia). 


El sistema estaba preparado para realizar inspecciones manuales y barridos. En el caso 
de estos últimos, por medio de una interfaz, se controlaban dos grados de libertad (X, 
Y). Los barridos automatizados podían ser tipo A-scan o B-scan y también tenía 
posibilidad de ajustar manualmente el eje Z donde se encontraba montado el palpador 
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Sistema SENDAS (4) 


Interfaz a Usuario 


camas po [ear] son] ro ca] Ala 
poa | bee (GmerarAB]sea |ro [ca |etcos | 
lis 
Retardo Gate 0 1 9 4 5 6 


3 


fanoo 4] 100 


Premuestras Posmuestras 80 
Disparo Nivel E E 

= 
cc | [E E 
Función Frecuencia 


PEF o A? 
Filtro jo H uN : 


FT Video AF G : 

anancia 0 
[Video MF , 10 92 195 297 400 502 
” Wideo BF hi 6 F dBs | Muestras 


al 


Sistema para END ULTRASCOPE (1) 


Desarrollado en el año 2004 por la empresa DASEL Sistemas — Madrid, España. Es un sistema 
basado en la arquitectura AMPLIA, desarrollada dentro del proyecto MAGIAS (Métodos 
Avanzados de Generación de Imágenes Acústicas) del CYTED y orientado a la generación, 
adquisición y procesamiento de señal ultrasónica para END. 

Una de sus características principales es la integración de múltiples funciones de procesamiento 
digital de señal, en una única FPGA XC25200 de Xilinx, con una alta tasa de utilización de sus 
recursos (96%), algunas de las cuales son exclusivas de este sistema. 


2 UltraScope 
e view Tools Configuratico 


'm +/0:/0,! ¿[83 
Men | Puser | Recebrer | 
Channel 10 == 
iv Pulser Enable 
Y A-Scan Enable 

Receive 
* Pulse-Echo 


” Transmission 


« fntemal 
C Extemal 
c Encoder 
Cc PRF 


Transducer Frequency 


Man | Gates | Encodes | TabSheet | Anale 
Range Gain 


> E y +0 IO 


Propagation Path 


Delay Material Velocity Reject Attenuator “ Depth (z) 
| 10.0us ] 0tus DNA > MEMBED > 2 6 | com. 58 
4) » J »] 41) »] 


Sistema para END ULTRASCOPE (2) 


Principales Funciones 


Función Característica 


Filtro FIR 31 coeficientes de 8 bits, programable. Opera en tiempo real (80 MS/s) 
Filtro EMI Elimina ruido impulsivo producido por interferencias EMI 

Curva CAT Compensación Atenuación-Tiempo: 512 puntos, programable, 48 dB 
Diezmado Programable 1:1 a 16:1 (equivalente a la frecuencia de muestreo) 

Compresión de Programable 1:1 a 128:1, sin pérdida de información (amplitud-tiempo) 

trazas 


Número de canales De l a 8 canales multiplexados (1 80 MS/s, pulso-eco o transmisión 


Detección de Rectificado doble onda y filtro IR, programable. 

envolvente 

Ventanas 3, pueden ser independientes o ligadas, operar con RF o Video 
MONOPICO 


Ventana multipico — 1, registra hasta 128 picos con umbral y ancho de validación programables 
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Esquema de un Sistema Analógico para A-scan 


ALMACENAMIENTO 
DE INFORMACION 


CONVERSOR 
AD 


RECEPTOR> 


PULSER 


CONVERTIDOR 
DIGITAL DE 
BARRIDO 


—<p—— PALPADOR 


TEJIDO 


jaa ORGANO 
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Relación entre A-scan y B-scan 
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Correspondencia entre A-scan y B-scan (Modo A y Modo B) 


Análisis de la relación entre A-scan y el B-scan obtenido a partir del primero. Es 
importante observar que la parte de mayor brillo en el registro del B-scan corresponde 
a valores de menor amplitud del A-scan. Inverso a la definición de B-scan. 

Se trata de un criterio de codificación de la información. El B-scan, por su 
naturaleza, se basa en el concepto de brillo : mayor amplitud, mayor brillo 


Modos 
Palpador 
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Vista de B-scan 
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Diferencias entre Modo A, B y M 


MODO E 


MODO M 
Ñ e 
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Correspondencia entre A-scan y B-scan y Modo M 


tiempo 
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B-scan: Probetas con taladros de fondo plano y horizontales 


a b | 
Material: Aluminio —= SS 2 
+ 


==... -----o-o- 


30 1] l l l e | 
50 100 150 200 250 300 350 400 


50 100 150 200 250 300 


B-Scan 


Alta resolución del B-Scan con solo dos datos por disparo 
¿Por qué la diferencia del fondo en ambas imágenes? 
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Control Atenuación Distancia — CAD (1) 


El control de ganancia visto, afecta a toda la traza, es decir, todo el eco es afectado 
por el valor de ganancia seleccionado. Existen algunas aplicaciones donde es 
necesario amplificar más un sector en particular de la traza y existen otras donde se 
necesita dar un perfil de ganancia diferente a distintos sectores de la traza. Los 
equipos modernos permiten al operador generar una curva de ganancia arbitraria, 
este control especial se lo denomina control atenuación distancia o CAD. que se 
genera con la asistencia de un mouse. 


Amplificador con Ganancia 
Controlada por Tensión 


Eco > o> Eco amplificado 


<Ú DAC K*= Ganancia Digital Generada 
por el mouse 


El número de puntos de la curva CAD puede estar en el orden de 128 o 256. 
En estos sistemas, la curva CAD generada es empleada para variar la ganancia del 
amplificador de entrada del receptor de ultrasonido 
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Control Atenuación Distancia - CAD (2) 


Ejemplo de Curva CAD compuesta por una ganancia fija de 10 dB mas una 
trapezoidal de 40 dB. Ecos procesados de un envase cerrado de puré de tomate 


con CAD triangular 
Envolvente A-scan 


A-scan final 


MN 
| 
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Control Atenuación Distancia — CAD (3) 


Ejemplo de Curva CAD compuesta por una ganancia fija de O dB mas una 
trapezoidal de 50 dB. Ecos procesados de un envase cerrado de puré de tomate 


con CAD triangular 


A-scan 


5v 


--5V 


5V 


o 


us 270 


Curva CAD 


Envolvente A-scan 


0 us 270 


Envolvente A-scan Final ei 


Otros Métodos Usados en END (1) 


1- Configuración de palpadores pitch-catch 
pitch-catch (lanzar y capturar). Este montaje de palpadores se usa entre otros, para la 
caracterización de defectos en núcleo de materiales compuestos. el ángulo de disparo del 
ultrasonido sea tal, que se generan ondas transversales. En emplean palpadores angulares y 
en general ondas de superficie y ondas reptantes o laterales (creeping waves) que viajan por la 
superficie y por el espesor de la pieza 


2- Técnica CHIME 
Basada en pitch-catch, permite trabajar con combinaciones de las ondas arriba citadas y 
detectar defectos basados en señales de difracción, da una idea del volumen donde se 
producen los efectos de difracción 


Para estos métodos el sistema de registro debe estar configurado en pulso — 
transparencia 
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Otros Métodos usados en END (2) 


3- Técnica TODF 

TOFD o Tiempo de transito por difracción (Time Of Fligh Difracction). Se basa en la 
medida del tiempo de transito de señales de generadas entre ambos palpadores por el 
defecto. 


Las señales de difracción están precedidas por una onda 
superficial y seguidas de la reflexión en la superficie 
opuesta a la reflexión en el fondo del material 
inspeccionado 


Eco de fondo 


Onda superficial 


Defecto, parte 
superior 


Defecto, parte 
inferior 
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Otros Métodos usados en END (3) 
Cálculo del tamaño y profundidad del defecto 
Secuencia de pasos 


1- Tiempo de transito de la onda lateral 
tt, == (1) 


2- Tiempo de transito sobre el borde superior del 
defecto debido a la difracción 


Eco de fondo 


JA D? + Ss? O) 5- Distancia a la que se encuentra el defecto 
AS E respecto de la superficie de inspección 
UA (5) 
3- Tiempo de transito sobre el borde inferior del D= 2 Ct 8 


defecto debido a la difracción 6- Altura del Defecto 


po 6) L= [06-57 -D (6) 


7- Espesor de la pieza 


4- Tiempo de transito debido a la reflexión en el 


1 (7) 
fondo de la pieza H ==. /¡Cf?, —S* 
$ 


4H*+S* (4) 
ly = 8- Altura del Defecto 


C 
s=,C?,-4P Y 2, 


Otros Métodos usados en END (4) 


Código de Simulación para el cálculo del tamaño y profundidad del defecto 


%todf 

%eEstima la altura L del defecto y la distancia D normal a la superficie de la pieza, en 
%función de la velocidad del sonido en el material, el tiempo transito t1 (parte superior del 
%defecto), t2 (parte inferior del defecto 

%tef (reflexión sobre el fondo de la pieza y separación de los palpadores. Además indica 
%el tiempo de transito de la onda lateral tl y el espesor H. 


%de la pieza. Los tiempos citados se obtienen del correspondiente A-scan 
function [tl, D,L,H]=tofd(t1,t2,tef,C,S) 

tl=S/C; 

D=0.5*sqrt(C12*t112-512); 

L=0.5*sqrt(C*2*t222-512)-D; 9 


H=0.5*sqrt(C12*tef"2-S12)-D; 


Ejemplo 7.2 

Aplicar el método TODF a una pieza de acero de 350 mm por 340 mm con un espesor de 13 
mm. Los palpadores de incidencia angular de 60 grados, separados a 96 mm. La c de las 
ondas longitudinales es de 5920 m/s. Del A-scan se tiene que el tiempo de transito de la onda 
lateral es de 16.7 us, el tiempo desde el borde superior del defecto es de 16.51 us, el del 
borde inferior de 17.25 us y el tiempo de transito debido a la reflexión en el fondo de la pieza 
es de 16.8 us. 

De acuerdo a los datos: [tl,D,L,H] = tofd(16.51e-6,17.25e-6,16.8e-6,5920,96e-3) 


resulta: tl = 1.62 us, D = 9.2 mm, L = 8.2 mm, H = 3.8 mm 
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Tema Ill 
Transductores de Ultrasonido, Lentes Acústicas y Arrays 
Primera parte 
Transductores de Ultrasonido y Lentes Acústicas 


TEMA 


Definición de transductor 

Modelo y Respuesta Mecánica de un Transductor para Ultrasonido 

Tipos de Sensores y sus Características Relevantes 

Modelo Eléctrico de un Transductor Piezoeléctrico 

El Campo de un Oscilador Ultrasónico 

Excitación Pulsada y Resolución Axial 

Concepto de Lentes y de Arrays Acústicos 

Características constructivas y variantes de Transductores Piezoeléctricos 

Características Eléctricas, Acústicas y Frecuenciales de Transductores 
Piezoeléctricos 
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Referencias 


Transductores de Ultrasonido y Lentes Acústicas (1) 


Definición 

En general el nombre de sensor, se refiere a dispositivos diseñados para cuantificar 
o detectar parámetros específicos, por medio de elementos transductores. En un 
sensor, los transductores son los elementos que desarrollan la función de 
transformación de una magnitud física en otra 


Emisión 


ELECTRICO >| MECANICO ACUSTICO 


Recepción 


ELECTRICO | MECANICO F | ACUSTICO 


79 


Transductores de Ultrasonido y Lentes Acústicas (2) 


Modelo y Respuesta Mecánica de un Transductor para Ultrasonido 


Ñ Rm 
M EE + in +Kx=0 | =-—— | Coeficiente de amortiguamiento 
dt dt 2M 


Tiempo de relajación ? = 


de)=x, e“ colw,t-0) + 


Transductores de Ultrasonido y Lentes Acústicas (3) 


Tipos de Sensores y sus Características Relevantes 
Existen diversas tecnologías para realizar esta función y en líneas generales 
podemos citar las siguientes: 


Electrostáticos — Electromecánicos - Piezoeléctricos 
Transductores capacitivos. 
Su empleo está limitado a aplicaciones en aire. 


ranuras 


aluminio 


celofán 


a b 
Frecuencia de resonancia 


R 1 
Coeficiente de amortiguamiento a-E-1 f= n T, 
| P 


E 


TP: tensión aplicada, |: ancho de las ranuras, p, densidad de la lámina, n: 
una constante que del modo de vibrar de la membrana 81 


Transductores de Ultrasonido y Lentes Acústicas (4) 


Transductores acústicos electromagnéticos (EMAT Electro Magnetic Acoustic 
Transducer) 


Aprovechan la ley de fuerzas F a las que se ve sometida una corriente eléctrica | en un 
campo magnético B tal como lo establece la ley de Lorentz 


Transductores Piezoeléctricos 
Los transductores piezoeléctricos se basan en el efecto que lleva su nombre 
descubierto en 1880 por J. y P. Curie. Cuando a un material piezoeléctrico se 
lo somete a un esfuerzo mecánico, se generan potenciales eléctricos entre 
sus caras. Estos mismos materiales muestran el fenómeno inverso, esto es, al 
aplicarles un potencial eléctrico sufren una deformación mecánica (ideales 
para pulso-eco). La frecuencia de oscilación se sitúa en el rango 5 kHz a 100 
MHZ. 
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Transductores de Ultrasonido y Lentes Acústicas (5) 


Modelo Eléctrico de un Transductor Piezoeléctrico 

Las cerámicas piezoeléctricas presentan una amplia variedad en los valores de las 
constantes elásticas, eléctricas y piezoeléctricas. Además, los materiales cerámicos se 
pueden fabricar en una extensa gama de formas y tamaños y su eje de polarización 
puede orientarse según la geometría y el modo de vibración deseado. Actualmente, los 
materiales cerámicos más conocidos y utilizados, son los que se conocen con el 
nombre de PZT constituidos por titanatos zirconatos de plomo. Existen otras 
composiciones comerciales basadas en titanatos de bario, metaniobatos de plomo, 
niobatos de sodio, etc. La posición lider de los PZT se debe a su potente efecto 
piezoeléctrico y elevado punto de Curie, junto a la amplia gama de propiedades que 
ofrecen con pequeños cambios en su composición 


dd 1 
=C a 1) AA NE 
(Es) ñ C, Qx f.,,) 


34 36 38 40 
Frecuencia [kHZ] 


_L2xf wL 
Q Q 83 


Transductores de Ultrasonido y Lentes Acústicas (6) 


El Campo de un Oscilador Ultrasónico 

Hasta aquí se ha supuesto que los frentes de onda ultrasónicos son uniformes. En la 
realidad varían por la geometría de la fuente, la distancia y orientación relativa entre el 
material y el oscilador, el medio de acoplamiento entre la fuente y el material, etc. El 
campo ultrasónico producido por un oscilador, es equivalente al que se produciría al 
atravesar un frente plano un orificio cuya geometría coincide con la del oscilador, y que 
está rodeado por una lámina infinita perfectamente rígida. Se puede demostrar que 
traspasada la lámina, se crea un frente plano limitado por el perímetro de la apertura 
que se propaga de forma uniforme, y, unas ondas, comúnmente llamadas ondas de 
borde, que surgen invertidas desde el perímetro de la apertura y se propagan como 
ondas libres en el espacio 


El campo creado por el oscilador, es el resultado de las 
interferencias entre el frente plano y las ondas de borde. 
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Transductores de Ultrasonido y Lentes Acústicas (7) 


sen TT D? 1 D = diámetro 
P=Bi———-—  delacerámica , 
SA  Z|  Z=distancia 


Z= 0/8) 2= 0/4). Z 
Vista espacial 


Campo cercano 


A 
a Apertura del haz  |y, = e [erados | 
A A ——— A 


i i 

| 0 D*/4), > Campo lejano 
j 

| 
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Transductores de Ultrasonido y Lentes Acústicas (8) 
Código de simulación 


function [Rel, N, fhi, P]=pre_emi_disco(D,lambda,n,Po) Ejemplo 

%Cálculo del campo cercano D= 24mm A=1.5 mm 
N=D*2/(4*lambda); z=.1:.1:n*N; cte=(D/2)12; 
%Cálculo de P a partir de la Ecu. (3.39) 


P=2*Pc*abs(sin((pi/laambda)*(sqrt(cte+z."2)-z))); 
%Cálculo de P simplificado, Ecu. (3.44) 
cte1=pi*D"2/(8*lambda); 
P1=Pcero*2*abs(sin(ctel./z)); 


1. — 
D/A=8, N=48. y =85" " ] 
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Transductores de Ultrasonido y Lentes Acústicas (9) 


Radiación de un Transductor tipo Pistón 


z, distancia primer cero J,(+3.83)=0 


J,(X), función de Bessel de 
primer orden 


ángulo de divergencia 
nOs 


=1 sobre el eje z z 
Y, =sen” (1215) 


Lóbulo principal Lóbulos secundarios 


0.9 


07 


POS 


Patrón de radiación Polar 87 


Transductores de Ultrasonido y Lentes Acústicas (10) 


Código de simulación 


function pre_rad_pcir1(D,lambda,gamma_min,gamma_max) 
%Conversión de grados a radianes 

fac=pi/180; 

gamma_min=gamma_min*fac; 
gamma_max=gamma_max*fac; 


rel=D/lambda; 
gamma=-gamma_min:.0017:g3amma_max; 
X=pi*rel*sin(gamma); 
P=abs(2*Besselj(1,X)./X); 


Ejemplo 
D/A = 16, D/A = 8 y D/A = 4. 


Línea continua D/A = 16, 
línea con puntos D/A = 8, 
línea con punto y raya D/A = 4. 


Angulo de apertura 
[Grados] 
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Transductores de Ultrasonido y Lentes Acústicas (11) 


Excitación Pulsada y Resolución Axial 


Señal de eco S)= Rf (1) mít) | p" -A 


e 
Portadora Rf(t) = sen [2 TT Jo f E 6] 7 
Parte ascendente > Función Potencial = tf" os PN 
. . - 0 
Parte descendente > Función Exponencial = e ** D 05 1 15 2 25 3 35 4 
A ——W 
tl 


Envolvente/Modulante mit) — A, a e - 


Modelo Polinómico del eco 
n = 2, para f, <= 320 KHz — n = 3, para f, > 320 KHz 89 


Transductores de Ultrasonido y Lentes Acústicas (12) 


La presión a la salida de un transductor es proporcional a la energía almacenada 
en sus oscilaciones. Se puede demostrar: 


Wo Surge una pregunta: ¿es deseable tener un 
E transductor con alto o bajo Q? 
pS En e Evidentemente depende de la aplicación. 


Si se trabaja con técnicas pulsadas, pulso-eco, conviene que el Q sea bajo. La 
conveniencia de tener un transductor de bajo Q está relacionada con la denominada 
Resolución Axial. Por resolución axial se entiende la capacidad de medida del 
instrumento usado para resolver o detectar dos objetos cercanos (poca 
separación entre ellos) sobre el eje axial o eje de la dirección en la que es emitida la 
onda de presión del transductor. 

% y 


/ Dirección del haz 


Dirección Axial 


E y 


Dirección Lateral 


Con otras palabras: Resolución axial es la capacidad para discriminar dos 
reflectores que se encuentran alineados con el eje del transductor 
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Transductores de Ultrasonido y Lentes Acústicas (13) 


Resolución axial o capacidad para discriminar dos reflectores que se 
encuentran alineados con el eje del transductor. Este parámetro depende 
del ancho de banda del pulso 


La siguiente figura muestra los ecos procedentes de las caras 1 y 2 de una 
placa estrecha. Los ecos resultan claramente distinguibles cuando se utilizan 
pulsos de corta duración, en cambio, con pulsos más largos, los ecos 
provenientes de las caras de la placa se superponen, por lo que se dificulta la 
evaluación 
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Transductores de Ultrasonido y Lentes Acústicas (14) 


Resolución Axial y Lateral 


Resolución Axial Resolución Lateral 


Depende del Q y de 1 Depende de D, del punto de enfoque y de 1 


Ejemplo 

Se dispone de un palpador que opera a una frecuencia de 3.5 MHz con un Q de 
7, un diámetro de 1.5 cm, la longitud focal del lente es de 5 cm y la longitud de 
onda en el medio es de 0.043 cm. Determinar la resolución axial y lateral 


_Q0,_1 _ 
= 4 A= di = 0.075cm Siel Q se duplica 


5 RA== 1=20.043= 0.15cm 
RL=2.44 joo: =0.35cm 
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Transductores de Ultrasonido y Lentes Acústicas (15) 


Resolución lateral o capacidad para separar dos puntos que están en una recta 
perpendicular a la dirección de propagación de los ultrasonidos. La resolución 
lateral viene dada por el ancho A del frente acústico utilizado. Para puntos del campo 
lejano el ancho del haz acústico aumenta linealmente con la distancia al transductor 


En el primer caso, la separación entre los puntos es inferior al ancho A del haz, por 
lo que resultan indistinguibles en el corte lateral. En cambio, cuando la distancia 
entre los puntos es superior al ancho del haz, éstos resultan perfectamente 
observables en la imagen. 
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Transductores de Ultrasonido y Lentes Acústicas (16) 


Las lentes acústicas realizan el mismo efecto de las lentes ópticas, es decir, 
modifican el perfil del haz ultrasónico con ánimo de obtener mejoras en los 
parámetros de las imágenes, especialmente la resolución lateral, intensidad de las 
señales ultrasónicas, uniformidad de las señales, etc. Existen numerosos perfiles de 
lentes, pero entre ellos destacaremos los dos tipos más usados: la lente esférica 
convergente y la lente cónica. 


Profundidad del Foco 


e Diámetro del Punto Focal 


Y 


ES | Foco 


eS 


Transductores de Ultrasonido y Lentes Acústicas (17) 


Ejemplo 

Se dispone de un palpador que opera a una frecuencia de 3.5 MHz con un Q de 7, un 
diámetro de 1.5 cm, la longitud focal del lente es de 5 cm y la longitud de onda en el 
medio es de 0.043 cm. Determinar la profundidad del foco, el diámetro del punto focal 
y relacionarlos con la resolución axial y lateral 


PF =3 (E 0.043=1.43cm 
A | 15 
' 0 ' 
¡ ds | d=2.44 (E 0.043=0.344cm 
pases A AR Es | 

l A —————>> 
1 RA ¡ 0.7 
| 1 RA=% 1=0.043=0.075cm 
A 4% 4 

ás RL= 244/0043 = 0.35 cm 
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Transductores de Ultrasonido y Lentes Acústicas (18) 


Modificación del lóbulo de directividad 


Transformador acústico 
Modifica el lóbulo de directividad en forma PASIVA basado en la teoría de Líneas 


de Transmisión Exponenciales 


op 1080p 10 p 
PEE 
dx" Soxox ví of” 


Arreglos acústicos (Array) 
Modifican el lóbulo de directividad en forma ACTIVA a 
través de RETARDOS emulando a una LENTE OPTICA 


Array lineal 


Array circular, Dcha: Array sectorial Array cuadrado bidimensional 
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Ejemplo de aplicación de ultrasonido en aire en el campo de la 
Asistencia en la movilidad de Ciegos 


Trabajo de Fin de Carrera — Ingenieria Electronica 
Facultad de Ingeniería — UNSJ - 1999 


Disposición de los Sensores 
Bastón Anteojos 


60” 


El control de la emisión y 
recepción de cada uno de 
los tres transductores que 
emulaban un array, estaba 

y cargo de un 
microcontrolador tipo PIC 


-60* 
-90* 


Este trabajo obtuvo el premio al Tecno-Emprendedor otorgado por el BBVA Banco 
Francés en el año 1999 y fue presentado en IBERSENSOR 2000, Madrid, España 98 


Implementación 


Características constructivas y variantes de Transductores 


Piezoeléctricos (1) 


1-Transductores para aire y medios gaseosos 


Estructura Interna PIE ASIOnES 


Malla metálica Resina epoxi 


O: 
NE E 


A e 
2-0.64 +0.1 =] 
S 
5.0+0.3 
MA40S4R/S MA405S5 
ho ul * : ElAJ Code 
O: Rors 
in mm 


Resistentes al agua 


MA40E7R/S 
*: ElAJ Code 
D:Rors 
in mm 


[o] 
WE 
Ss 
bx] 


y 
m Sa 
(9) YA 
y [ NN vi Es 


AS 
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Características constructivas y variantes de Transductores 


Piezoeléctricos (2) 
2-Transductores para END y Medicina 


2-1- Incidencia Normal 
Focalizados y No Focalizados 2-2- Incidencia Angular 


Cerámica 


y Amortiguador 
Protección 


Cerámica 


V540-SM 


Lente ABWM-5T-X? 


+ 


Protección 


C540-SM 
ABSA-5T-X? 
dl Protección 


2-3- Para Doppler continuo 


Cerámica 


Protección | 101 


Características Eléctricas, Acústicas y Frecuenciales de 


Transductores Piezoeléctricos (3) 
Características Eléctricas y Frecuenciales 


1- Impedancia Eléctrica en Función de la Frecuencia 


MA40B8R MA40S4R 


40 


8 
Sensitivity (dB) 
g € 


30 35 40 45 50 30 35 40 
Frequency (kHz) Frequency (kHz) 


2- Patrón de Radiación 


MA40B8R MA40S4R 
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Características Eléctricas, Acústicas y Frecuenciales de 


Transductores Piezoeléctricos (3) 


2-4- Características Temporales 


2-5- Características Frecuenciales 
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Aplicaciones Industriales del Ultrasonido 
Ultrasonido para Uso Médico 


Prof. Dr. Angel C. Veca 
Marzo 2018 


Tema lll 
Segunda parte 
Arrays Ultrasónicos y Focalización Electrónica 


TEMA 


Un poco de Historia 

Concepto de Array 

Arrays lineales de N fuentes equiespaciadas 

Sistemas de imagen ultrasónica 

Focalización y deflexión del haz 

Técnicas electrónicas de focalización dinámica en recepción 

Elementos de simulación 

Referencias 106 


Un poco de Historia 


El concepto de array nace en el ámbito de la aviación alrededor de 1927, como 
guía de aterrizaje en condiciones de baja o nula visibilidad . A partir de un 
conjunto de 2 antenas se formaba un lóbulo de radiación llamado trébol de 4 
hojas. La aeronave cuya ruta coincidía con el espacio entre las hojas del trébol 
tenia asegurado el descenso seguro 


Evolución de los Barridos Mecánicos usados en Diagnóstico Médico que llevaron a la 


necesidad del array 
Codificador 
Codificador de posición 


de Posición 
7 
y” pam Sistema de Motor 
ya engranajes 
E ! h Palpado: 
Ati 1 SS r 
stático nao ! 
Operador Pda | 
: y 1 y s Palpador 
1 N ; ¡ A 
Palpador ———»- / A a 
4 k 4 - 17 s 
| s / I Ús 
od / Oscilante j , 
Material bajo estudio ———=»> ; y : 
», l eS / i z 
A] te A 4 ze 


! “  Rotatorio  * 
En ninguno de los casos se modificaba el diagrama de Ñ | á 


directividad del palpador, se variaba mecánicamente SÓlO — "rematan 
la posición XY j 107 


Concepto de Array 


Un array es un conjunto de transductores idénticos ordenados según una cierta geometría y 
alimentados para obtener un diagrama de radiación predefinido. En este contexto la palabra array 
se utilizará para designar a un arreglo de cualquier geometría. Sin embargo, los fundamentos 
físicos se explicarán en base a un array lineal formado por varios elementos piezoeléctricos 
dispuestos a lo largo de una línea y separados regularmente por una distancia finita 


La expresión de la onda de presión generada 

por una de las fuentes en un punto del Cuando se 

espacio está dada por trata de dos 
fuentes 


Arrays lineales de N fuentes equiespaciadas 


A 

/1l Lóbulo principal 
¡¿Lóbulo de rejilla 
A 


| Lóbulos secundarios / | 
| | X 


H(8): función de directividad o factor de array 


LL 
nivel de los lóbulos laterales (NLL) NLLldB]|= SES d*N: apertura del array 


Lóbulos de Rejilla 
LR 


nivel de los lóbulos de Rejilla (NLL) |NLR|dB]|= 201085 
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Parámetros Clave para la Calidad de la Imagen 


RESOLUCIÓN AXIAL ANCHO DE BANDA DISPLAY 
RESOLUCIÓN LATERAL TAMAÑO APERTURA FOCALIZACIÓN 
RANGO DINÁMICO, RSR LÓBULOS SECUNDARIOS RESOLUCIÓN TEMPORAL 


OTROS 
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Concepto de Focalización Electrónica (1) 
Focalización: Mejora la resolución lateral y la RSR 


MONOELEMENTO 


FOCALIZADO cacaos 


RETARDOS 
AÑADIDOS 


RETRASOS 
PROPAGACIÓN 


Q 


ELEMENTOS 
VIRTUALES 


Y 


IDEALMENTE: Focalizar en cada punto de la imagen, en EMISIÓN y RECEPCIÓN 
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Concepto de Focalización Electrónica (2) 


Focalización en Emisión 


RETARDOS 


| RETARDOS 
DIGITALES PULSERS 


FRENTE DE ONDAS 


ELEMENTOS 


VIRTUALES Ñ 


alla 
Dbboor 


| 
LL 


PROGRAM. PROGRAM. 
AARARO RETARDOS AMPLITUD 
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Concepto de Focalización Electrónica (3) 


Focalización en Recepción 


No existe limitación física en cuanto al número de focos por línea (como en emisión) 


De hecho, se puede “mover” el foco siguiendo la propagación del ultrasonido: 


ARRAY 
N e Sentido de propagación 
dE A A A A  _ —_——> 


Con un único disparo se puede 
focalizar en TODOS los puntos 


E Pe Llegada ecos F1 Ñ 
N 5 de la línea... 


Llegada ecos F2 


El único problema es electrónico: ¿Cómo variar DINÁMICAMENTE los retardos? 
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Focalización sobre el Eje Principal 


TIEMPO DE PROPAGACIÓN (Tp): 


1 2 
Í, Ll 


C C E Z 


APROXIMACIÓN BINOMIAL (x << z) 


2 


TIEMPO DE PROPAGACION DESDE EL 
ELEMENTO CENTRAL AL FOCO 


DIFERENCIA EN TIEMPO DE PROPAGACIÓN 
DESDE EL ELEMENTO ACTUAL (x,0) ALFOCO 


NON ANNA 
N 


T, no depende 
de x 


donde T, = constante suficientemente grande para que T, > 0. 
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Focalización y Deflexión del Haz 


=> — A 
Ez Ea x"  xsenO 
z a z a r=xR“+x?*-—2xRsen0 = R|1+ ES 
2R?—OR 
a 
E [] 
[] 
[] 
"A A 
[8 
x[m / 0 
M >, 
D z 
O Z DEFLEXION 
[] 
[E 
O x xsenO 
INSTANTES DE DISPARO PARA M7” —7T — E 
- FOCALIZAR Y DEFLECTAR X A De A 
[] 
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Ejemplo 


ARRAY de N=64 elementos, d=1/2=0*25 mm (D=16 mm), FOCO a R=830 mm. 


c= 1500 m/s (medios biológicos), deflexión Oy ,x = 45” 


Los retardos se calculan en tiempo no real y se programan en la electrónica. 
TT x ñ xsenO 0 A 
* "2% 2Re C E A A 


T, (ns) 0=02 258 
T, (ns) 0=452 


Gran Rango de retardos (varios Us ) 


Elevada Resolución (pocos ns) 
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¿Qué Significa Focalización Dinámica? 


1. Una operación sencilla: Sumar el valor de N muestras para cada foco 


2. Las muestras deben corresponder al “instante de llegada” desde cada foco a 
cada uno de los elementos del array. Para ello, la señal recibida por cada 
elemento se retrasa una cantidad variable con el tiempo, función de x, R, 0 


3. Las señales deben ser de RF o contener información de fase para producir 
interferencias constructivas (donde interesa) y destructivas (donde no), lo que 
demanda una alta precisión en los retardos (una fracción del periodo de RF) 


4. La secuencia de salida (resultado de la suma) es un A-scan en RF, del que hay 
que extraer la envolvente (intensidad) y preparar para su visualización 


5. Generalmente este proceso debe realizarse en tiempo real 
Ejemplo: imagen sectorial con: 


128 líneas, R = 1024 muestras/línea, array lineal N = 64, a 
| = 25 imágenes/s 
Capacidad de procesamiento (1)+(2)+(3): 


L*R*N*I = 128*1024*64*25 = 210 MS/s 
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Focalización Dinámica en Emisión y Recepción 


Si se dispone de focalización dinámica en recepción y capacidad de 
deflexión electrónica también se puede aproximar a la focalización 
dinámica en emisión usando las Técnicas Confocales 


Focalizan en emisión y recepción, operando en tiempo real: 
* Múltiples focos/línea en emisión y recepción 


* Tiempo equivalente a 1 disparo/línea 


3 disparos casi simultáneos Aprovecha la profundidad de campo del foco 
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Técnicas Confocales, Principio de Funcionamiento 


Rotación de los disparos múltiples DISPARO 1 


= 
pao) 
N 


Haz 


En consecuencia: Con K disparos se obtienen K zonas focales 
que cubren K líneas 119 


Resumen de Focalización 


1. Los arrays permiten focalizar en EMISIÓN y en RECEPCIÓN utilizando lentes 
arbitrarias. 


2. La focalización en EMISIÓN es un proceso relativamente sencillo (retrasar una 
señal digital). 


3. La focalización en RECEPCIÓN es un proceso mucho más complejo (retrasar 
múltiples señales analógicas, conservando sus características: forma, amplitud, 
frecuencias, etc.). 


4. En general se requieren grandes aperturas para mejorar la resolución lateral 
(sen a. =1/D) 


5. Grandes aperturas = pocos elementos muy distanciados (económico) = 
LOBULOS DE REJILLA muchos elementos a d = 1/2 (costoso y complejo) 
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Sistemas Electrónicos de imagen ultrasónica 


Dos aproximaciones 
Barrido lineal Barrido sectorial 


121 


122 


Sistemas de imagen ultrasónica Sectorial 
(phased arrays) 
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Sistemas de imagen ultrasónica Sectorial 


(phased arrays) 


Realizan la deflexión del haz 


Muchos elementos activos (todos) 
Adquisición en coordenadas polares 


Cubren mayor área con menor apertura 


Su mayor inconveniente es la complejidad de la electrónica 
que requieren 
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Sistemas de barrido lineal (analógicos) 


N ELEMENTOS, M<N CANALES 


xM 
IMAGEN 
CONTROL 
ga 
ARRA 

DETECC. E 

ENVOLV. 
xM 


Sistemas muy simples 


Ejemplo: N=128 elementos 
M= 8 (pulsers, amplificadores) 
1 único A/D. 
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Sistemas de barrido lineal con Focalización Dinámica 


N ELEMENTOS, M<N CANALES 


CONTROL 


FOCALIZACIÓN 
- EN RECEPCIÓN DETECC. 
(POCOS ENVOLV. 
CANALES) 


Ejemplo: N=128 elementos 


M= $ (pulsers, amplificadores) 


M=8 A/D+ FOC.REC. 
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Sistemas de barrido sectorial (phased arrays) 


N ELEMENTOS, N CANALES 


IMAGEN 
CONTROL 
ARRA 

A/D 

DETECC. SCAN 
A/D ENVOLV. CONV. 
A/D 

Ejemplo: N=64 elementos 


N= 64 (pulsers, amplificadores) 
N= 64 A/D+ FOC.REC. 


Sistemas complejos 127 


Algunas cifras... 


a ii 


: RECEPCIÓN | 
| EMISIÓN — |! 


TRANSDUCTOR ¡ ADQUISICIÓN | ¡FOCALIZACIÓN: 


PANTALLA 


Array de N 
elementos 


N = 32 a >500 Decenas de MHz Varias GS/s 


Tiempo real sin artefactos 


Resolución axial y temporal 


Resolución lateral sin lóbulos de rejilla 
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Realidad 


Sólo para conseguir una alta resolución de retardos (rango dinámico): 


Ejemplo: f, = 5 MHz, fmu = 80 MHz !! ... unas $ veces por encima de Nyquist 


Si 200 us/línea (15 cm en medios biológicos) en cada canal se adquieren: 


M= 200-80 = 16.000 muestras/canal (para cada línea) 


MEMORIA datos CONVERSOR A/D FOCALIZADOR SISTEMA 


2x1MGO 12ns Alto coste Alta velocidad 


Después de todo este esfuerzo, sólo se utiliza una pequeña fracción de las muestras 
adquiridas: aquellas que se suman para cada foco (1024 en cada canal) 


MÁS DEL 90% DE LA INFORMACIÓN ADQUIRIDA ES INÚTIL 


CONCLUSIÓN: hay que utilizar métodos alternativos más eficientes 


Objetivo: Reducir la frecuencia de muestreo manteniendo una alta resolución temporal 
en los retardos 


Problemas a resolver: 

1- Reducir la cantidad de memoria de focalización 

2- Reducir la cantidad de memoria de señal adquirida 

3- Reducir la velocidad de muestreo, manteniendo una alta resolución temporal 
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Técnicas electrónicas de focalización dinámica en recepción 


1. 
Ze 
de 
4, 
3: 


Técnicas analógicas 

Técnicas de Apertura Sintética (SAFT) 
Técnicas digitales convencionales 
Muestreo paso-banda e interpolación 


Técnicas de Muestreo Selectivo (SST) 
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Técnicas Analógicas de Retardo Variable 


1. MEDIANTE LÍNEAS DE RETARDO CONMUTABLES: 


SES 


2. MEDIANTE LÍNEAS DE RETARDO VARIABLES (CCDs): 


LEO TAO AT 
E E TT E LL LE J 


EIA == MÍA == ARAN 
E E E LL LL LE _JL 


3. OTRAS TÉCNICAS (HETERODINACIÓN HACIA ARRIBA Y USO DE FILTROS SAW) 
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Técnicas Mixtas 


1. LÍNEAS DE RETARDO => PEQUEÑOS RETARDOS (<1 us) 
2. RETARDOS DIGITALES: PROGRAMABLES. 
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Técnica de Apertura Sintética (SAFT) 


N disparos y adquisiciones 
para una imagen 


PULSER 


Composición de señales por 
software o por hardware 


IT IT III 


133 


Focalización Básica por Software 


A Crear una tabla de índices k a muestras correspondientes de cada canal y sumar. 


profundidad z (para una línea) A — VA 


aaa A A 
ES AF 3 
5] [261] A d, 

ÓN EN 


| ll 
JJ MEMORIA DE 
FOCALIZACIÓN 


¿Cuánta memoria hace falta F 
para a los índices? ; Ñ 
N ÓN R puntos de AE a Í : vi 
en N-R índices por cada línea. Ejemplo: ¿Cuánta memoria hace falta 
A 2 
N=64, R=1024, L=128 > 2% punteros =83M para almacenar las señales* 


. 1 
con un tiempo de acceso = 4 ns Ñ 
O a+? dilo 


Fr.  ————> 
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Fundamentos de la Focalización Digital en Tiempo Real 


A Tamer] Memoria de señal 
> LU 
a CARTESIANAS 


(RASTER) 
POLARES 


A-SCAN (RF) 
y DET. CONV. 
EM DISPLAY 


LA RESOLUCIÓN DEL RETARDO ES 
IGUALAL PERIODO DE MUESTREO: 


ALTA FRECUENCIA DE MUESTREO. 
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Focalización con Interpolación 


ZA GRUESO (periodo de muestreo verificando Nyquist) 


CAJD Laos IN 
HE 1 ep FINO 


»nQn in 
ES BA EA 


Xx 
- A 
Fa d INTERP y ds 
CONV. 
E dd SCAN 
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Repaso de las Bases Teóricas para la Interpolación 


Sea x(k) = x(kT) la secuencia que resulta al muestrear una señal continua x(t) con 
periodo de muestreo T que verifica el criterio de Nyquist: T<1/(2fmax). 


x(t) Ñ Me x(k) 


Según la teoría de muestreo, puede reconstruirse EXACTAMENTE x(t) a partir de x(k)... 


iria La reconstrucción por interpolación 
LINEAL no es muy “EXACTA”... 


sinc(Y/LAT/T) | xp(k) 


== AT 
x(t- AT) 
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La interpolación no es trivial... 


ZA GRUESO (periodo de muestreo verificando Nyquist) 


CAJD COIN, 
HE 1 e FINO 


»nQn 
ES BA EA 


+ 
XK Ñ 
Fa O INTERP ) 


Ej. Interpolación de 3 puntos: 


CiXa + Cox + CiXu1 


Los coeficientes C,, C, y C, 
varían dinámicamente 


137 


Técnicas de Muestreo Selectivo (TMS) 


En lugar de RETRASAR las señales, se basan en adquirir las muestras en el 
instante de llegada desde cada foco a cada elemento, almacenarlas y sumarlas 
ordenadamente. 

El RELOJ de muestreo de cada canal, de frecuencia NO UNIFORME, determina 
los instantes de adquisición: 


Z 
X 


2 
Tozg lx, R,0)=T, O 
Rc  cC 
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Focalización Encadenada 


Profundidad del FIFO 
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Focalización Encadenada: Operación 


LINEA DE FOCOS 
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Arquitectura TMS para Focalización 


REALIZACIÓN DISTRIBUIDA DE FIFOs 


Hi 1 2 4 8 


yz, 
2 
2 


<||S<||S||S||S| ||| 
O||U||O|[O||U||U||U 


OU 


ANNAN 


> 
pa 


Ejemplo: para un array con N=64=2%8 elementos, se requieren 384 posiciones en vez 
de 4096 
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JN QA Un bb Y 


VYVVVYVVYVVYVWVY 


Focalizador TMS en AMPLIA (proyecto MAGIAS) 


Módulo anterior Módulo siguiente 
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Introducción a la Simulación de Arrays (1) 


%C= velocidad del sonido en el medio considerado, f0 = frecuencia del emisor, fm = 
%frecuencia de muestreo, k = número de onda, gN = número de puntos de discretización del 
espacio 

function [p]=monopolo(np,C,fo,fm) 


% grilla de puntos 
x0=0; 

y0=np/2; 
[x,y]=meshgrid(-np/2:np/2); 

jwt=sqrt(-1)* 2*pi*f* 0; 

k = 2*pi*f/c; 

r = sqrt(x.22+y.12); 

p1 = exp(jwt-k*r); 

%Graficacion 

imagesc(abs(p1)), axis equal; title(Monopolo') 


Monopolo 


monopolo(400,344,100); ma), 
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Introducción a la Simulación de Arrays (2) 


Dos fuentes monopolares 


function [x,y,p1,p2] = bipolo2(c,fo,Lx,np,n,sep) 
deltaphi=n*pi/2; % diferencia de fase entre p1 y p2 
lambda = c/fo; % longitud de onda 
x0=sep*lambda; 


Fm = np/Lx; % frecuencia de muestreo del espacio 
[x,y]=meshgrid(xmin:1/Fm:xmax); % grilla de puntos 

wt= 2*pi*fo*0; 

k = 2*pi*f/c; 

r = sgqrt((x+x0/2).22+(y+y0/2).22); 

p1 = exp(1j*(wt- k*r + 0*pi))./r; % monopolo 

rp = sqrt((x-x0/2).22+(y-y0/2).22); 

p2 = exp(1j*(wt-k*rp + deltaphi))./rp; % monopolo desplazado 
p = p1+p2; 

figure(1); imagesc(xmin:1/Fm:xmax,xmin:1/Fm:xmax,log(abs(p))) 
axis equal; title('Bipolo'); xlabel('X”); ylabel('Y”); 
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Introducción a la Simulación de Arrays (3) 


Dos fuentes monopolares 


Simular el campo de presión para el caso de dos fuentes monopolares separadas a una de 0.34, 0.54, 0.64, y A 
que oscilan a una frecuencia de 40 kHz, sumergidas en aire y diferencia de fase igual a cero. 


0,14 
bipolo2(344,40€e3,0.2,2003,0,0.1); 


A E Pol 
bipolo2(344,40€e3,0.2,2003,0,0.5); -0.06 
bipolo2(344,40€3,0.2,2003,0,1); 0.04 
Debido a que se supone que el medioes 
aire y teniendo en que cuenta la e 
frecuencia de trabajo es de 40 kHz, se 0.08. 


adopto Lx en 20 cm 


Introducción a la Simulación de Arrays (4) 


Focalización 
function [phi ] =array(N,sepy,R,c,f,Lx,condir,gN,dy,xe) 
xe =0; 
| = 3e-3; % ancho del transductor (m) 
lambda = c/f; % longitud de onda 
deltay=sepy*lambda; 
xmin = -Lx/4; xmax = Lx*3/4; ymin = -Lx/2; ymax = Lx/2; 
Fse = gN/Lx; % frecuencia de muestreo del espacio 
[x,y]=meshgrid(xmin:1/Fse:xmax,ymin:1/Fse:ymax); % grilla de puntos 
w= 2*pi*f; k = 2*pi*f/c; 
p=x*0; % inicializo la presión en cada punto 
for ¡=1:N 
r = sqrt((x-xe).22+(y+deltay*N/2-deltay*i).22); % distancia cada punto al elemento ¡ 
thet = atan2((y+deltay*N/2-deltay*i),(x-xe)); 
If condir 

dir = ( sin(pi*llambda*sin(thet))./(pi*llambda*sin(thet)) ); 
else 

dir=1; 
end 
phi(i) = (sqrt(R(1).22+(R(2)-deltay*i+deltay*N/2).22)-sqrt(R(1).224+R(2).22))/c; 
p = p + dir.*exp(1j*(w*phi(i) - k*r ))./r; % acumulador 
end 
figure(1) 
imagesc(xmin:1/Fse:xmax,ymin:1/Fse:ymax,log(abs(p))) 
axis equal; hold on 
plot(R(1),R(2),'or); plot(0,(-N/2:N/2)*deltay,'.b') 
hold off 
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Introducción a la Simulación de Arrays (5) 


Simular la distribución del campo acústico de un array de 4 elementos apoyado sobre una placa 
de aluminio, que trabaja a una frecuencia de 1 MHz. a- Focalizando sobre el eje. b- Focalizando 
en el punto R=[0.1 0,05]. c- Para las mismas condiciones, para 16 elementos focalizando sobre el 
eje y d- Focalizando para el punto R=[0.1 0,05]. 


array(4,[0.1 0.0]); 

array(4,[0.1 0.05]); 

array(16,[0.0 0.0]); 
( 


array(16,[0.0 0.05]); 
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Aplicaciones Industriales del Ultrasonido 
Ultrasonido para Uso Médico 


Prof. Dr. Angel C. Veca 
Marzo 2018 


Tema IV 
Simulación de la Señal de Eco 


TEMARIO 


Modelación y Simulación de Trazas A-scan usando Matlab 
Parámetros Característicos de la Envolvente de un Eco 
Generación de una Traza Simple 

Generación de una Traza Múltiple 

Detección de Máximos de una Traza 


Referencias 
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Modelación de trazas (1) 


Motivación 
A partir del modelo de eco con envolvente polinómica: 
REO) =C, (t—T,)"e270) sen[27x0, (+-T,)+ 0] ute) 
1- Generar una traza de un eco teniendo en cuenta: 
Parámetros del palpador, Frecuencia de muestreo 


Tiempo de transito, Relación Señal/Ruido 
Amplificación 


2- Generar una traza de tres ecos teniendo en cuenta: 
los parámetros anteriores 


3- Aplicar a una traza generada un algoritmo de multipicos 


Se usa el modelo polinómico debido a que la mayoría de 
los palpadores comerciales responden a el (51 


Modelación de trazas (2) 


Determinación del Tiempo de Transito 
A-sscan Rf(t)=C, (t-11) "2%" sen [2iro, (t—1t)+ p] utt) 


Envolvente de Rí(t) mt) = C, (1 o tt) n p7li=n) 


Donde: Co es la amplitud máxima, ala constante de amortiguamiento del 
palpador, n constante del modelo polinómico con valor = 2 para frecuencias 
menores o iguales a 320KHz y 3 para frecuencias mayores. 

tt es el tiempo de transito, pla fase del eco, wO la frecuencia angular de 
trabajo y u(t) la función escalón unitario 
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Modelación de trazas (3) 


Calculando la derivada primera de Rí(t) respecto de (t-tt) e igualando a cero 


Derivando por segunda vez respecto del tiempo la señal de eco e igualando 
a cero, se obtienen los puntos de inflexión: 


mit) 


A vn Punto de inflexión Ascendente 


194 
= 458 vn Punto de inflexión Descendente 
(Y Mis 


ia 


Mid 


Lg 
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Generación de una traza simple con Matlab (1) 


function [eco,env]=genera_poly15(f0_nor,t,alfa,nm,EX,snr) 


donde: eco es el A-scan, env, envolvente del A-scan 
f0_nor, frecuencia de muestreo normalizada 

t, vector que define la longitud del registro 

alfa, constante de amortiguamiento del palpador 

nm, número de muestras al máximo del eco 

EX, amplitud del eco 

snr, relación señal ruido en dB 


f0_nor, frecuencia de muestreo normalizada: 


Relación entre la frecuencia de trabajo del palpador fo, respecto a la 
de muestreo fm 
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Generación de una traza simple con Matlab (2) 


Cómo se relaciona el número de muestras con el tiempo? 


l 
tiempo Tm(s) = l A 


nm= : = 
Frecuencia de muestro fm 


Periodo de muestro 


Estimando la distancia a recorrer por el ultrasonido, las características del 
medio y la frecuencia de muestreo, puede determinarse el número de 
muestras al máximo de la envolvente y el número de muestras totales de la 
simulación. El siguiente ejemplo aclara lo citado. 


Ejemplo 

Se dispone de una pieza de aluminio (c = 6000 m/s) de 2 cm de espesor 
para ser inspeccionada usando un palpador de 5 MHz, con constante de 
amortiguamiento de 0.13 y se muestrea a una frecuencia de 66 MHZ. 


Calcular: 
El tt teórico 
El nm al máximo de la envolvente 


La longitud máxima del registro . 155 


Generación de una traza simple con Matlab (3) 


tt =2* espesor/C =2*0.02/6000 = 6.66us 


Tm=1/ fm=1/66-10* =15.1ns 
nm, = tt/Tm= 6.66us/15.8ns = 440 muestras 


t.. =6.66:-10%+3.49-107 =7us 
nm... =£,.,. /Tm= Tus/15.8ns = 463 muestras 


Para visualizar la señal de eco completa se puede adoptar un tiempo de 
4 veces la diferencia de tiempo t,,, ,-tt, debido a que se trata de una 
función exponencial negativa, con valor próximo a cero a partir de 4 


constantes de tiempo, por lo tanto, 


nm, = Ut, —14)/Tm= 4(Tus —6.66us)/15.8ns = 90muestras 


El número total de muestras del registro es de 553, quedando en la 
función genera_poly15, el vector de muestra en [0:552] 
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Generación de una traza simple con Matlab (4) 


genera poly15.m 


function [eco,env,to]=genera_poly15(f0_nor,t,a,n,nm,EX,snr,gra) 
% Modelo de envolvente polinómica 
w0O_nor=2*pi*f0_nor; 
tmax=round(n/a); 
b=(a*exp(1)/n).Mn; 
% Generación del eco 
tO=0; 
tx=[zeros(1,(t0+nm)-tmax) 1:length(t)-((t0+nm)-tmax)]; 
e=b*tx.m.*exp(-a*tx); 
rí=sin(wO_nor*tx); 
eco=EX*e.*rf; 
% Ruido aditivo 
y=10*(snr/20); 
rui=(rand(size(tx))-.5)/y; 
eco=eco+rul; 
%Extracción de la envolvente 
env=abs(hilbert(eco)); 
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Generación de una traza simple con Matlab (5) 


fm = 66 MHz, fo = 5MHz,t =[0:512], nm = 100, Ex = 1, SNR = 40aB, alfa = 
0.13 
[eco,env]=genera_poly15(5/66,[0:512],0.13,100,1,40,1) 


A-scan 


| | ] 
0 100 200 300 400 500 600 
Número de muestras 


Envolvente 


| 
0 100 200 300 400 500 600 
Número de muestras 
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Generación de una traza simple con Matlab (6) 
Cambiando el nm al máximo a 300 


[eco.envi=senera nolv15(5/66.10:5121.0.13.300.1.40.D 


4 


0.5 


A-scan 
o 


-0.5 


0 100 200 300 400 500 600 
Número de muestras 


Envolvente 


Ds ye A A A RS a | en AA Egel 
0 100 200 300 400 500 600 
Número de muestras 
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Generación de una traza simple con Matlab (7) 
Cambiando la snr de 50 a 20 dB 


[eco,env]=genera_poly15(5/66,[0:5 12],0.13,300,1,20,1) 


A-scan 


0 100 200 300 400 500 600 
Número de muestras 


Envolvente 


0 100 200 300 400 500 600 
Número de muestras 
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Generación de una traza simple con Matlab (8) 


Cambiando fm de 66 a 33 MHz 
[eco,env]=genera_poly15(5/33,[0:512],0.13,300,1,50,1) 


A-scan 
o 


| | | 
0 100 200 300 400 500 600 


Cambiando fm de 33 a 15 MHz 
[eco,env]=genera_poly15(5/ 15,[0:512],0.13,300,1,50,1) 


0.5F 


o Ale 


A-scan 


-0.5; 


4 


| | | | 
0 100 200 300 400 500 600 
Número de muestras 
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Generación de una traza de 3 ecos en Matlab (1) 


traza.m 


function [traza,env] = genera_traza(f0_nor,t,a,snr,EX1,EX2,EX3,nm1,nm2,nm3) 


donde: eco, A-scan, env, envolvente del A-scan 
f0_nor, frecuencia de muestreo normalizada 

t, vector que define la longitud del registro 

alfa, constante de amortiguamiento del palpador 


nmi, número de muestras al máximo del eco 


EX, amplitud del eco 
snr, relación señal ruido en dB 


Elementos comunes a la traza(f0_nor,t,a,snr,) 


Elementos particulares de cada eco(,EX1,EX2,EX3,nm1,nm2,nm3) 
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Generación de una traza de 3 ecos en Matlab (2) 


genera_traza.m 


function [traza,env] = 
genera_traza(f0_nor,t,a,n,snr,EX1,EX2,EX3,nm1,nm2,nm3,g) 
%Generacion de los ecos 
[eco1,env1]=genera_poly15(f0_nor,t,a,n,nm1,EX1,snr,0); 
[eco2,env2]=genera_poly15(f0_nor,t,a,n,nm2,EX2,snr,0); 
[eco3,env3]=genera_poly15(f0_nor,t,a,n,nm3,EX3,snr,0); 
% Generacion de la traza 

traza=eco1+ec02+ec03; 

env=env1+env2+env3; 

%Graficacion 

If g == 

figure(1);subplot(2,1,1);plot(traza) 
subplot(2,1,2):plot(env) 

end 
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Generación de una traza de 3 ecos en Matlab (3) 
snr = 50 dB 
genera_traza(f0_nor,talfa,5 0 100,200,300,1,0.7,.4); 


Amplitud 


Amplitud 


Traza de tres ecos 
4 


0.5 
0 
-0.5 
-1 
0 100 200 300 400 500 600 
Numero de muestras 
Envolvente de la traza de tres ecos 
1.5 
1 IM | 
0.5 
0 


Numero de muestras 


| Ss 
0 100 200 300 400 500 600 
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Generación de una traza de 3 ecos en Matlab (4) 


Cambiando nm3 de 300 a 400 muestras 
genera_traza(f0_nor,t,alfa,40, 100,200,400, 1,0.4,.7); 


Traza de tres ecos 
1 


0.5- 


0 


Amplitud 


-0.5- 


| 
0 100 200 300 400 500 600 
Numero de muestras 


Envolvente de la traza de tres ecos 


Amplitud 


0 100 200 300 400 500 600 
Numero de muestras 
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Generación de una traza de 3 ecos en Matlab (5) 


snr = 20 dB 


genera_traza(f0_nor,t,alfa,20,100,200,400,1,0.4,.7); 


Amplitud 


Amplitud 


4 


Traza de tres ecos 
T 


T 
| 
| 
0.5- Z 
0 a) 
-0.5- al 
-1 | | | 
0 100 200 300 400 500 600 
Numero de muestras 
Envolvente de la traza de tres ecos 
1.5 
1 lo. al 
0.5 
0 | | | 
0 100 200 300 400 500 600 


Numero de muestras 
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Detección de Múltiples Máximos con Matlab (1) 
function [C] = multipicos2015(A,1,M,u,p) 


Detector de múltiples picos de la traza A 

les un vector con los márgenes del algoritmo 
M es el intervalo del máximo 

u es el valor del umbral 

prueba, control de gráficos 


Multipicos a partir de un umbral 


0.5 


Amplitud 
o 
o 


. Amplitud 


0 50 100 150 200 250 300 350 le | | 
5 2 95 3 
Numero de muestras " nia MERA 130 zen 2. sa. 


Numero de muestras 
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Detección de Múltiples Máximos con Matlab (2) 
multipicos2015.m 


function [C] = multipicos2015(A,i,M,u,g) 
a =1(1);]=1; 
while (a < 1(2)) 

fA(a) > ([A(a-M:a-1) A(a+1:a+M)]) 
A(a); 
aja=a+M; ¡=j+!l; 


end 
end 
for i=1:max(size(B)) 

1f B(1,2)>u, B(1,2)=B(1,2); 

else 

B(1,2)=0; 

end 
end 
E=find(B(:,2)>u); C=[B(E,1) B(E,2)); 
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Detección de Múltiples Máximos con Matlab (3) 


Amplitud del eco a detectar en la tra<a, >= 0,88 


multipicos2015(eco,[90 450],10,0.88,1); 


Traza de tres ecos 


0.5 
Jo) 
? 
a 0 y 
3] 
<X 
0.5 | 
-1 | 
0 100 200 300 400 500 600 
Numero de muestras 
4 
0.5- 
0 [AAN 
-0.5 
a] | 
0 100 200 300 400 500 600 
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Detección de Múltiples Máximos con Matlab (4) 


Amplitud del eco a detectar en la tra<a, >= 0,85 


multipicos2015(eco,[90 450],10,0.85,1); 


Traza de tres ecos 
T 


Amplitud 
[e) 


| 
0 100 200 300 400 500 600 
Numero de muestras 


600 
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Detección de Múltiples Máximos con Matlab (5) 


Amplitud del eco a detectar en la tra<a, >= 0,5 


multipicos2015(eco,[90 450],10,0.5,1); 


Traza de tres ecos 


0.5- 
Jo) 
2 
a ] 
É 0 
<X 

0.5- Z] 

1 | 

0 100 200 300 400 500 600 
Numero de muestras 
j 0 100 200 300 400 500 600 
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Tema VI 
Efecto Doppler en Medicina 


TEMARIO 


Concepto de Efecto Doppler Mapas de Flujo 

Concepto de Doppler Direccional y No Direccional Comparación entre Doppler Color y Doppler 
Modulación y Demodulación de AM de Potencia 

Doppler Continuo Ejemplos de combinación Doppler-Imagen 
Demodulación Coherente (No Direccional) Referencias 


Demodulación Heterodina (Direccional) 

Proceso de Visualización 

Concepto de Doppler Pulsado 

Operación del Doppler Pulsado 

Resumen de Diferencias entre Doppler Continuo 

y Doppler Pulsado 

Algunos Ejemplos 

Frecuencia de Muestreo 

Limitaciones del Doppler Pulsado 195 


Concepto de Efecto Doppler (1) 


|] 
Sistema de medición ' Sistema de traducción 
|] 


Universo del sistema 
S 


Velocidad 


Frecuencia 


y 
Sentido 
del objeto 


Un > Us 
V>f>A 
V>oA 


donde: V está representando la velocidad del objeto incluyendo el sentido del 
movimiento (directo e inverso), f el cambio de frecuencia y Á la amplitud proporcional 
a la velocidad 


196 


Concepto de Efecto Doppler (2) 


Como se vio en el tema l, de los tres casos posibles, interesa sólo el caso 3: Emisor en 
reposo y receptor moviéndose en la misma dirección y cualquier sentido 


Ñ pa 
f E a fa 
El ) ) ) ) [ro Ecuación general del Efecto Doppler sn [H4)s 


c+V, 


ve 
Para este caso V¿ = 0 Fr E 2) E Operando fi. = R Í, E dl E 
C 


Definiendo la frecuencia diferencia o frecuencia Doppler, como fy = fa — fe Je le Í, 
En pulso eco LS 2, 
- E 
Haciendo el análisis para el caso de la ... 2 fecos O 


hemodinámica, se arriba a Sa y C 
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Concepto de Efecto Doppler (3) 


EJEMPLO 1: 
Frec. emisión f¿ = 1 MHz 
Velocidad Va = 1 m/s 
c = 1500 m/s 
fp = 1,3 KHz 


EJEMPLO 3: 
Frec. emisión f¿ = 5 MHz 
Velocidad Va = 1 m/s 
Cc = 1500 m/s 
fp = 6,7 KHz 


EJEMPLO 2: 
Frec. emisión f¿ = 2 MHz 
Velocidad V¿= 1 m/s 
c = 1500 m/s 
fp = 2,66 KHz 


EJEMPLO 4: 
Frec. emisión f¿ = 10 MHz 
Velocidad Va = 1 m/s 
c = 1500 m/s 
fp = 13,33 KHz 
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Ayudas al Diagnóstico 


$4 


ESTENOSIS 


FLUJO 
DIRECTO 
FLUJO 
INVERSO 


SEÑAL DOPPLER 


BIDIRECCIONAL 
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Revision de Modulación y Demodulación de AM 


Modulación FO o F(w) FO) 


cos w,t F(0)cos m1 5[F(w+ w, )+F(w=w, )] 


ty — aldo => ((t) cos w.! 


f(t)cos”w.to FO) + os 2w,)+ Flw-2w,)| 


F(w) 


Demodulación 


cos wi ene | F(1)+ f()cos2w 1] 


f(t) cos w,t 2 
Ñ ¡ y Multiplicador HA COS Filtro Pasa 10) 
X Bajos 
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Concepto de Doppler Direccional y NO Direccional 


Sistema 
Receptor 
(Traducción) 


CW Direccional 


CW NO Direccional 


Visualización Direccional 


Z(f) VELOCIDAD 


(FRECUENCIA) 


3 KHz 8 KHz 


AMPLITUD 


cm/s 


VELOCIDAD SANGUÍNEA 
EN EL TIEMPO 


A $ E 
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Doppler Continuo CW 


EMISOR 
QA 
Ss Ss 
N= XX 
Pp? 


R¿cos(w,t+wWst+0f) + Racos(w,t-Wat+0p) 
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Demodulación Coherencte (no direccional) 


R(D=Rocos(wt+dp) + Recos(wotewpt+9p) + Racos(Wwytatró) 


OSCILADOR 


FILTRO P.BAJO 


200 Hz 5 KHz 
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Demodulación Heterodina (Direccional) 


OSCILADOR 


3 KHz 8 KHz 


FILTRO P.BAJO 
(2'9 KHz) 


FILTRO P.ALTO 
A) 


3 KHz 8 KHz 


3 KHz 8 KHz 
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Decodificación en Ccuadratura 


OSCILADOR 


1/2 


D(t) 


Q(t) 


D(t) = A cos(0,) + Bf Cos(Wt + 07) + Bp COS(Wpt - 0p) 
Q(t) = A sen(0)) + B, sen(w,t + ds) - Bah sen(wat - dp) 
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Extracción de las Componentes del Flujo (1) 


D(t) = A cos(0,) + Bf Cos(Wt + 07) + Bp COS(Wpt - 0p) 


Q(t) = A sen(4)) + B, sen(wst + ds) - Bah sen(wat - 0) 


C.C. FORWARD REVERSE 


D(t) 
Bf sen(wpt - dp) 


B; cos(wft + 0f) 
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Extracción de las Componentes del Flujo (2) 


Técnica Digital 


D(t) = B, cos(wft + df) + Br COS(Wat - da) = Xp(t) + Xp(t) 
Q(t) = B, sen(wft + df) - Bg sen(wrt - 0a) = H[X£(t)] - H[Xp(t)] 


D'(t) = H[Xp(t)] + HIXp(t)] , mn» D'(t) + Q(t) = 2 H[X;(t)] 
Q'(t) = HIXf(t)] - H[Xp(t)] D'(t) - Q'(t) = 2 H[Xp(t)] 


D(t) HILBERT 
(FIR orden N) 2Bp sen(wat - 


2B, sen (wpt + dp) 
RETARDO 


ZN 
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Extracción de las Componentes del Flujo (3) 


Uso de Tablas 


D(t) = B, cos(wft + df) + Bp COS(Wat - dp) 
Q(t) = B sen(wft + df) - Ba Sen(wgt - (A) 
Para cada instante t: 
| X(t) = f(D, Q) = B, cos(wrt + d+) 


Y (t) = g (D, Q) = Ba sen(wat - dp) 


D(t) 
TABLA 


f(D,Q) pH)» X(t) 


NN 


Q(t) 


> VO 
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Extracción de las Componentes del Flujo (4) 


Uso de la Transformada de Fourier 


Este método se basa en las propiedades de simetría de la transformada compleja de Fourier 
donde, si la señal de entrada está en cuadratura, las componentes directas e inversas del 
flujo aparecen a ambos lados del lóbulo principal de la transformada 


D 0.1 (az D.3 D.4 0.5 0.65 o 
Tiempo (s) 


Resultado del uso de la Transformada de Fourier donde se 
aprecia una mala elección de la frecuencia de muestreo 
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Proceso de Visualización 


SEÑAL DOPPLER 


RELOJ MUESTREO 


Z(f) 


da CONVERSOR ele 
: 


3 KHz 8 KHz 
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Concepto de Doppler Pulsado - PW 


RANGO 


DETECCION DE FASE 


de /dt 


ZI2 


1/PRE —> 


FFT 


Operación del PW 
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PW No direccional: Esquema en Bloques y Operación (1) 


La señal a la salida del oscilador maestro (OM) pasa a través de una llave al driver de 
salida, que es controlada por la CPU. El intervalo de tiempo que la llave permanece 
cerrada para que pase un número completo de ciclos del OM depende del muestreo 
del sistema. El pulso resultante de RF es introducido a un amplificador de radio 
frecuencia y enviado al transductor emisor. Los ecos que se producen en el tejido 
retornan y son captados por el transductor receptor, son amplificados (RF) y se 
introducen en el mezclador X para realizar la demodulación. A continuación pasan por 
un filtro pasa bajos que elimina la componente de RF y son muestreados con la llave 
L) la que se encuentra abierta durante la transmisión. Esta llave es en realidad un 
circuito de muestreo y retención. A la salida del muestreador hay un filtro anti aliasing 
(pasa banda), para suprimir las señales espurias que introduce el muestreo. 
Finalmente, un amplificador de baja frecuencia da los niveles de ganancia necesarios. 
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PW No direccional: Esquema en Bloques y Operación (2) 


Ventanas de tiempo en PW El diagrama muestra el pulso de transmisión de 
duración t,, la ventana de tiempo tg en la que se 
encuentra trabajando la llave L, y la ventana 
temporal desde que se habilita la transmisión y se 
abre L,. Definiendo con d, la distancia 
comprendida entre el transductor y el inicio del 
rango, entonces, 


cl La -1,) 
d, = 7 
2 
clt,+t,) 
donde: c es la velocidad del sonido en el medio. De forma similar se define d,  d, = > 


como la distancia entre el transductor al final del rango de medición 


Por lo tanto el rango puede ser expresado como: d, = d,— d, reemplazado, se tiene: 


_ Cl, +, 


" 2 


donde: d, es la demora entre la transmisión y la habilitación de L,. Se puede demostrar que si 
se tiene la mejor sensibilidad para una resolución axial una dada. Para t, > t, la sensibilidad 
decrece debido a la posibilidad de cambios producidos por scatering de señales con oposición 
de fase 
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Ejemplo PW No direccional 


Para un sistema PW que trabaja a 5 MHZ a una profundidad inicial de 1 cm en un medio con 
c = 1600 m/s. Determinar los tiempos de trabajo definidos en las Ventanas de tiempo en PW 
(Transparencia anterior). 


Como f, = 5 MHz, el periodo t, = 0,2 us. Suponiendo que se adopta que el pulso de 
excitación tenga un ancho t, = 10 t,, t, = 2 us. Teniendo en cuenta lo visto, si se toma t, = t,, 
ty = 2 Us. 

Por otra parte, el tiempo necesario para penetrar 1 cm, es de 6.25 us, entonces t¿ = 6,25 us. 
Finalmente, para el periodo de repetición del ciclo, PRF se tiene en cuenta que, después de 


la primera reflexión si se espera un tiempo igual a t,, no existirán ecos por la atenuación del 
medio por lo que PRF será: 


PRF =21(t, +21,)=29us 
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PW No direccional: Estimación de la velocidad y el ángulo de 


incidencia 


Si se graficq en función de la frecuencia de trabajo y 
de la velocidad del sonido en el medio en cuestión, 
obtenemos para 1, 2 y 5 MHz las siguientes gráficas: 


5 MHz 2 MHz 


Además de la frecuencia de trabajo, el ángulo con que el 3 
haz de ultrasonido incide sobre el tejido influirá en f.. 
Evidentemente v cos O será máximo para 8 = 0* y cero 

para O = 90* tal como lo indica la Fig. 10.16. El [ms] 
transductor no puede colocarse a 8 = 0%, entonces, debe 
buscarse el menor ángulo posible para obtener la mayor dd 


sensibilidad. 0<9<20" 
0<0<20" "| ma 


El cómputo de la velocidad se realiza con la expresión ya o. 2 80  Blgrados] 


vista 
772 ESC 
2 f cosO 


Si mantenemos constante c, f,, cos € y definimos una constante como cte = c/2 f, cos 6, 


entonces, 
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PW direccional: Esquema en Bloques y Operación (1) 


Lave 


> a IP Demodulador > B 


Básicamente se trata del sistema no direccional al que se ha agregado un demodulador, una 
unidad de muestreo y retención y un filtro pasa banda (antialiasing). Las unidades de 
muestreo y retención reemplazan a la llave L, . Los filtros de hasta 8 polos están 
encuadrados en una ventana frecuencial comprendida entre 100 y 800 Hz. Estos sistemas 
posen gran bW e alto slew rate. 


Las señales provenientes de cuerpos en estado estacionario y de pequeñas señales de los 
mismos en movimiento son amplificadas con la misma ganancia, el receptor requiere de un 
gran rango dinámico para evitar saturación de señales provenientes de las paredes de las 
arterias. Por ello, el receptor debe tener también la capacidad de recuperarse rápidamente de 
la saturación y disponer de un ancho de banda que es inversamente proporcional a la 
resolución axial del sistema. 
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Limitaciones del PW Direccional (1) 


El Doppler pulsado presenta dos grandes limitaciones: en primer lugar, igual que el CW la 
desviación Doppler es proporcional al coseno del ángulo que forma el haz ultrasónico con el 
flujo sanguíneo. Si el haz es paralelo al flujo, el ángulo es cero y el coseno 1, de manera que 
la desviación Doppler representa verdaderamente la velocidad del flujo. Con otros ángulos, 
esta premisa ya no es cierta. En la práctica no se utilizan ángulos mayores de 20* ya que 
producen estimaciones del flujo significativamente más bajos de lo que son en realidad. 


Ejemplo 

Estimar el error cometido en una medición cuando el ángulo que forma el haz ultrasónico y 
el flujo es de 2, 5, 20 30 y 45 grados. El instrumento está calibrado para la condición ideal, 
8 = O grados. Para esta condición, V = 0.5002 m/s, f¿ = 3.126 kHz, f, = 5 MHz y c = 1.600 
m/s. Por lo visto, se tiene: 


V,, =0,5005, V,, =0,5021, V,, =0,5323, V,. =0,5775, V,., =0,707 


En función de estas cifras, se deja que el alumno estime el error cometido para 6 > 0. 

En segundo lugar, la máxima velocidad del flujo que puede medirse viene determinada por la 
frecuencia con la que el sistema emite un pulso y abre la ventana de medición (t, en la Fig. y 
por la profundidad a la que se encuentra el volumen de muestra. A frecuencias bajas, el 
volumen de muestra es pequeño, por lo que, la posibilidad de encontrar anormalidades en la 
velocidad y dirección de los eritrocitos es escasa. A frecuencias mayores puede ocurrir que un 
segundo impulso sea emitido antes de que el eco del primero haya llegado de retorno al 
transductor. Por otra parte, cuando mayor sea la profundidad a la que se encuentra la 
muestra, más tiempo tarda el eco en volver, por lo que puede ocurrir que, la señal de retorno 
sea interferida por un nuevo pulso antes de que pueda ser captada por el transductor. 
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Limitaciones del PW Direccional (2) 


En la práctica, la máxima velocidad del flujo sanguíneo que puede ser determinada sin 
ambigúedad, viene dada por el límite de Nyquist que está dado por la frecuencia de repetición 
de pulso y la profundidad de la muestra. Si este límite se supera, se produce una saturación de 
la señal o "aliasing". El aliasing ocurre a partir de los 0.40-0.60 m/s [4,6,7], siendo así que el 
flujo cardiaco normal puede alcanzar 1.5 m/seg y en algunas patologías hasta 6 m/s. En estos 
casos el PW no tiene utilidad y hay que recurrir al CW 


Ambigúedad en la velocidad y en el rango 
Desplazamiento Doppler mad comparar la relación de fase con la señal de referencia OM. 
El máximo cambio de fase que puede obtenerse entre dos pulsos está limitado al intervalo -TT a 


+TT radianes (las medidas angulares se repiten cada 211 radianes). Por lo tanto la frecuencia de 
muestreo debe ser por lo menos el doble de la componente de frecuencia más alta. 


donde: fama, es el máximo desplazamiento Doppler que puede ser detectado 
Pomo = da sin ambigúedad y f, la frecuencia de repetición del pulso o frecuencia de 
EN muestreo. V..., sin ambiguedad, es: 
ds E fo = F. 2f Va, cosO 
d max A O O 


2 C 


La posición de la llave L, está determinada por el tiempo de demora, entre la transmisión y el 
comienzo del tiempo de adquisición de la señal recibida. Hay ambigúedad en la determinación 
del rango debido a que las señales que arriban al transductor, pueden ser ecos resultantes del 
pulso de transmisión anterior, localizándose alrededor de: cl ad ] 
dn = == 
2 


donde: n= 0 o entero positivo y ts tiene que ver con PRF 20 


Limitaciones del PW Direccional (3) 


Haciendo los reemplazos correspondientes, Amax Y Vmax a detectar sin ambigúedad: 


se trata de un sistema 


C 
d max == acoplado, donde se debe 
2 > ms 
: Se ti a lograr una situación de 
cf e tiene que — ' Y 
e s 8 f.cose| compromiso entre Vmax 
4 f,cosO y la profundidad máxima 
de trabajo. 
Vmax 
[mm/s] 
6 
En la práctica, los fabricantes de 
equipos los diseñan de acuerdo a 4 Le 
gráficas similares a ma), 
Es decir, Vmax = f(d) con f, como > NUS 
parámetro 5 
, 4 8 12 Distancia 
Ejemplo [mm] 


A partir de los datos obtenidos en el ejemplo anterior, estimar la d,,.. Y V,, para que no 
exista ambigúedad si f, = 5 MHz. Remplazando valores: 


fama = 17,24KHZz 
max =23,2cm 


Vox = 2,8m/s 2 


Resumen de las diferencias entre CW y PW 


l  Doppler Continuo (CW) | Doppler Pulsado (PW) 


Brinda distribución de velocidades de | Brinda distribución de velocidades 
Ventajas los glóbulos rojos de los glóbulos rojos 
Permite determinar velocidades Permite determinar velocidades 
máximas, mínimas promedio máximas, mínimas y promedio 
Presenta variación temporal del 
flujo 
Mejores posibilidades de acceso a El volumen de la muestra puede 
zonas pequeñas definirse a diferentes 
profundidades 
Permite estudiar sitios 


Emite y recibe continuamente 
determinados: resolución de rango 
Puede medir altas velocidades 9 


No permite estudiar sitios determinados: 


ambigúedad de rango ecocardiográfica para ubicar su 


posición 
Emite y espera el regreso del eco 
para hacerlo nuevamente 
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La velocidad máxima está limitada 
Desventajas por la profundidad: se produce 
aliasinc 
Siempre requiere información 


Algunos Ejemplos 


A la izquierda un eco doppler normal, a la derecha, un eco doppler anormal. Se 
observa, por un lado, el ensanchamiento de la base del registro y por otro, el 
incremento del reflujo 


Ejemplos de Eco doppler 


A A LL A A Md A 
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+” 


PLA IE. ARAFAT 


uu 
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Otras Limitaciones del PW Direccional 


La máxima velocidad del flujo que puede medirse viene determinada por la frecuencia 
con la que el sistema emite un pulso y abre la ventana de medición (td) y por la 
profundidad a la que se encuentra el volumen de muestra. 


A frecuencias bajas, el volumen de muestra es pequeño por lo que la posibilidad de 
encontrar anormalidades en la velocidad y dirección de los eritrocitos es escasa. 


A frecuencias mayores puede ocurrir que un segundo impulso sea emitido antes de 
que el eco del primero haya llegado de retorno al transductor. 


Por otra parte, cuando mayor sea la profundidad a la que se encuentra la muestra, 
más tiempo tarda el eco en volver, por lo que puede ocurrir que la señal de retorno sea 
interferida por un nuevo pulso antes de que pueda ser captada por el transductor. 


En la práctica, la máxima velocidad del flujo sanguíneo que puede ser determinada sin 
ambigúedad viene dada por el límite de Nyquist determinado por la frecuencia de 
repetición de pulso y la profundidad de la muestra. Si este límite se supera se produce 
una saturación de la señal o "aliasing". 
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Variantes de PW (1) 


Existen tres variantes básicas: 
* Doppler Color 

* Doppler de Energía 

* Doppler de Pared 


El Doppler Color (DC) utiliza el efecto Doppler 
generado por los glóbulos rojos en movimiento, 
utiliza el color para codificar la frecuencia de 
corrimiento Doppler. El color y los diferentes tonos 
de los mismos son determinados por la dirección y 
la magnitud del cambio de frecuencia. AMPLITUD 
Para codificar la turbulencia se usa la varianza, la 
que expresa las desviaciones de las diferentes 
velocidades presentes con respecto a la media. Y VARIANZA 
Como generalmente la profundidad y la resolución 
de velocidad son los parámetros más importantes, 
éstos se mantienen sacrificando el ángulo X 
sectorial, la frecuencia de imagen y la densidad de 

líneas. Los mapas en color del flujo se superponen 

sobre los ecocardiogramas que se muestran en 

gris. 


POTENCIA 
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Variantes de PW (2) 


El Doppler de Potencia no tiene problemas de aliasing ni dependencia angular, 
permite determinar el flujo de pequeños vasos y de baja velocidad. Recibe 
diferentes nombres: Color Power Doppler, Ultrasound Angio, Color Doppler Energy, 
Color Power Angio. Asigna valores a diferentes amplitudes de la señal Doppler. La 
amplitud está determinada por la concentración de los dispersores (eritrocitos) que 
producen el desplazamiento de frecuencia. La amplitud está relacionada con el 
volumen en movimiento. Integra el área bajo el espectro. La ventaja respecto a 
otros modos es la presentación uniforme del flujo, aún en aquellas zonas de 
incidencia normal. Sin embargo se pierde información de dirección, velocidad y 
característica del flujo. 


Distribución de Potencia Pseudo-Instantánea 


Potencia 


D.3 
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Frecuencia (Hz) Tiempo (s) 


Variantes de PW (3) 


Brinda distribución 
Ventajas de las velocidades espacial de las partículas 
en movimiento 
Permite calcular para cada 
ubicación el valor medio, la flujo 
varianza y el signo de la 
velocidad ángulo 


Desventajas No realiza una No permite evaluar 
caracterización completa | ninguna característica del 
del flujo flujo 


El Doppler Pared realiza el análisis del movimiento de las paredes cardíacas, 
donde se tiene baja velocidad y alta amplitud del flujo, a diferencia del resto del 
organismo donde se tiene alta velocidad y baja amplitud. Debido a que se trata de 
flujo turbulento y rotacional, el sistema dispone de filtros de amplitud para suprimir 
las señales producidas por el flujo, en general se tratan de filtros adaptivos. 
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Ejemplos de Combinación Doppler-Imagen 


PRF = 1.5 KHz 


- 19.2 
ecmís 
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Tema VI 
Efecto Doppler en Medicina 


TEMARIO 


Concepto de Efecto Doppler Mapas de Flujo 

Concepto de Doppler Direccional y No Direccional Comparación entre Doppler Color y Doppler 
Modulación y Demodulación de AM de Potencia 

Doppler Continuo Ejemplos de combinación Doppler-Imagen 
Demodulación Coherente (No Direccional) Referencias 


Demodulación Heterodina (Direccional) 

Proceso de Visualización 

Concepto de Doppler Pulsado 

Operación del Doppler Pulsado 

Resumen de Diferencias entre Doppler Continuo 

y Doppler Pulsado 

Algunos Ejemplos 

Frecuencia de Muestreo 

Limitaciones del Doppler Pulsado 195 


Concepto de Efecto Doppler (1) 


|] 
Sistema de medición ' Sistema de traducción 
|] 


Universo del sistema 
S 


Velocidad 


Frecuencia 


y 
Sentido 
del objeto 


Un > Us 
V>f>A 
V>oA 


donde: V está representando la velocidad del objeto incluyendo el sentido del 
movimiento (directo e inverso), f el cambio de frecuencia y Á la amplitud proporcional 
a la velocidad 
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Concepto de Efecto Doppler (2) 


Como se vio en el tema l, de los tres casos posibles, interesa sólo el caso 3: Emisor en 
reposo y receptor moviéndose en la misma dirección y cualquier sentido 


Ñ pa 
f E a fa 
El ) ) ) ) [ro Ecuación general del Efecto Doppler sn [H4)s 


c+V, 


ve 
Para este caso V¿ = 0 Fr E 2) E Operando fi. = R Í, E dl E 
C 


Definiendo la frecuencia diferencia o frecuencia Doppler, como fy = fa — fe Je le Í, 
En pulso eco ds Ry, 
- E 
Haciendo el análisis para el caso de la ... 2 fecos O 


hemodinámica, se arriba a Sa Ñ C 
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Concepto de Efecto Doppler (3) 


EJEMPLO 1: 
Frec. emisión f¿ = 1 MHz 
Velocidad Va = 1 m/s 
c = 1500 m/s 
fp = 1,3 KHz 


EJEMPLO 3: 
Frec. emisión f¿ = 5 MHz 
Velocidad Va = 1 m/s 
Cc = 1500 m/s 
fp = 6,7 KHz 


EJEMPLO 2: 
Frec. emisión f¿ = 2 MHz 
Velocidad V¿= 1 m/s 
c = 1500 m/s 
fp = 2,66 KHz 


EJEMPLO 4: 
Frec. emisión f¿ = 10 MHz 
Velocidad Va = 1 m/s 
c = 1500 m/s 
fp = 13,33 KHz 
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Ayudas al Diagnóstico 


$4 


ESTENOSIS 


FLUJO 
DIRECTO 
FLUJO 
INVERSO 


SEÑAL DOPPLER 


BIDIRECCIONAL 
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Revision de Modulación y Demodulación de AM 


Modulación FO o F(w) FO) 


cos w,t F(0)cos m1 5[F(w+ w, )+F(w=w, )] 


ty — aldo => ((t) cos w.! 


f(t)cos”w.to FO) + os 2w,)+ Flw-2w,)| 


F(w) 


Demodulación 


cos wi ene | F(1)+ f()cos2w 1] 


f(t) cos w,t 2 
Ñ ¡ y Multiplicador HA COS Filtro Pasa 10) 
X Bajos 
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Concepto de Doppler Direccional y NO Direccional 


Sistema 
Receptor 
(Traducción) 


CW Direccional 


CW NO Direccional 


Visualización Direccional 


Z(f) VELOCIDAD 


(FRECUENCIA) 


3 KHz 8 KHz 


AMPLITUD 


cm/s 


VELOCIDAD SANGUÍNEA 
EN EL TIEMPO 


A $ E 
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Doppler Continuo CW 


EMISOR 
QA 
Ss Ss 
N= XX 
Pp? 


R¿cos(w,t+wWst+0f) + Racos(w,t-Wat+0p) 
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Demodulación Coherencte (no direccional) 


R(D=Rocos(wt+dp) + Recos(wotewpt+9p) + Racos(Wwytatró) 


OSCILADOR 


FILTRO P.BAJO 


200 Hz 5 KHz 
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Demodulación Heterodina (Direccional) 


OSCILADOR 


3 KHz 8 KHz 


FILTRO P.BAJO 
(2'9 KHz) 


FILTRO P.ALTO 
A) 


3 KHz 8 KHz 


3 KHz 8 KHz 
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Decodificación en Ccuadratura 


OSCILADOR 


1/2 


D(t) 


Q(t) 


D(t) = A cos(0,) + Bf Cos(Wt + 07) + Bp COS(Wpt - 0p) 
Q(t) = A sen(0)) + B, sen(w,t + ds) - Bah sen(wat - dp) 
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Extracción de las Componentes del Flujo (1) 


D(t) = A cos(0,) + Bf Cos(Wt + 07) + Bp COS(Wpt - 0p) 


Q(t) = A sen(4)) + B, sen(wst + ds) - Bah sen(wat - 0) 


C.C. FORWARD REVERSE 


D(t) 
Bf sen(wpt - dp) 


B; cos(wft + 0f) 
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Extracción de las Componentes del Flujo (2) 


Técnica Digital 


D(t) = B, cos(wft + df) + Br COS(Wat - da) = Xp(t) + Xp(t) 
Q(t) = B, sen(wft + df) - Bg sen(wrt - 0a) = H[X£(t)] - H[Xp(t)] 


D'(t) = H[Xp(t)] + HIXp(t)] , mn» D'(t) + Q(t) = 2 H[X;(t)] 
Q'(t) = HIXf(t)] - H[Xp(t)] D'(t) - Q'(t) = 2 H[Xp(t)] 


D(t) HILBERT 
(FIR orden N) 2Bp sen(wat - 


2B, sen (wpt + dp) 
RETARDO 


ZN 
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Extracción de las Componentes del Flujo (3) 


Uso de Tablas 


D(t) = B, cos(wft + df) + Bp COS(Wat - dp) 
Q(t) = B sen(wft + df) - Ba Sen(wgt - (A) 
Para cada instante t: 
| X(t) = f(D, Q) = B, cos(wrt + d+) 


Y (t) = g (D, Q) = Ba sen(wat - dp) 


D(t) 
TABLA 


f(D,Q) pH)» X(t) 


NN 


Q(t) 


> VO 
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Extracción de las Componentes del Flujo (4) 


Uso de la Transformada de Fourier 


Este método se basa en las propiedades de simetría de la transformada compleja de Fourier 
donde, si la señal de entrada está en cuadratura, las componentes directas e inversas del 
flujo aparecen a ambos lados del lóbulo principal de la transformada 


D 0.1 (az D.3 D.4 0.5 0.65 o 
Tiempo (s) 


Resultado del uso de la Transformada de Fourier donde se 
aprecia una mala elección de la frecuencia de muestreo 
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Proceso de Visualización 


SEÑAL DOPPLER 


RELOJ MUESTREO 


Z(f) 


da CONVERSOR ele 
: 


3 KHz 8 KHz 
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Concepto de Doppler Pulsado - PW 


RANGO 


DETECCION DE FASE 


de /dt 


ZI2 


1/PRE —> 


FFT 


Operación del PW 
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PW No direccional: Esquema en Bloques y Operación (1) 


La señal a la salida del oscilador maestro (OM) pasa a través de una llave al driver de 
salida, que es controlada por la CPU. El intervalo de tiempo que la llave permanece 
cerrada para que pase un número completo de ciclos del OM depende del muestreo 
del sistema. El pulso resultante de RF es introducido a un amplificador de radio 
frecuencia y enviado al transductor emisor. Los ecos que se producen en el tejido 
retornan y son captados por el transductor receptor, son amplificados (RF) y se 
introducen en el mezclador X para realizar la demodulación. A continuación pasan por 
un filtro pasa bajos que elimina la componente de RF y son muestreados con la llave 
L) la que se encuentra abierta durante la transmisión. Esta llave es en realidad un 
circuito de muestreo y retención. A la salida del muestreador hay un filtro anti aliasing 
(pasa banda), para suprimir las señales espurias que introduce el muestreo. 
Finalmente, un amplificador de baja frecuencia da los niveles de ganancia necesarios. 
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PW No direccional: Esquema en Bloques y Operación (2) 


Ventanas de tiempo en PW El diagrama muestra el pulso de transmisión de 
duración t,, la ventana de tiempo tg en la que se 
encuentra trabajando la llave L, y la ventana 
temporal desde que se habilita la transmisión y se 
abre L,. Definiendo con d, la distancia 
comprendida entre el transductor y el inicio del 
rango, entonces, 


cl La -1,) 
d, = 7 
2 
clt,+t,) 
donde: c es la velocidad del sonido en el medio. De forma similar se define d,  d, = > 


como la distancia entre el transductor al final del rango de medición 


Por lo tanto el rango puede ser expresado como: d, = d,— d, reemplazado, se tiene: 


_ Cl, +, 


" 2 


donde: d, es la demora entre la transmisión y la habilitación de L,. Se puede demostrar que si 
se tiene la mejor sensibilidad para una resolución axial una dada. Para t, > t, la sensibilidad 
decrece debido a la posibilidad de cambios producidos por scatering de señales con oposición 
de fase 
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Ejemplo PW No direccional 


Para un sistema PW que trabaja a 5 MHZ a una profundidad inicial de 1 cm en un medio con 
c = 1600 m/s. Determinar los tiempos de trabajo definidos en las Ventanas de tiempo en PW 
(Transparencia anterior). 


Como f, = 5 MHz, el periodo t, = 0,2 us. Suponiendo que se adopta que el pulso de 
excitación tenga un ancho t, = 10 t,, t, = 2 us. Teniendo en cuenta lo visto, si se toma t, = t,, 
ty = 2 Us. 

Por otra parte, el tiempo necesario para penetrar 1 cm, es de 6.25 us, entonces t¿ = 6,25 us. 
Finalmente, para el periodo de repetición del ciclo, PRF se tiene en cuenta que, después de 


la primera reflexión si se espera un tiempo igual a t,, no existirán ecos por la atenuación del 
medio por lo que PRF será: 


PRF =21(t, +21,)=29us 
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PW No direccional: Estimación de la velocidad y el ángulo de 


incidencia 


Si se graficq en función de la frecuencia de trabajo y 
de la velocidad del sonido en el medio en cuestión, 
obtenemos para 1, 2 y 5 MHz las siguientes gráficas: 


5 MHz 2 MHz 


Además de la frecuencia de trabajo, el ángulo con que el 3 
haz de ultrasonido incide sobre el tejido influirá en f.. 
Evidentemente v cos O será máximo para 8 = 0* y cero 

para O = 90* tal como lo indica la Fig. 10.16. El [ms] 
transductor no puede colocarse a 8 = 0%, entonces, debe 
buscarse el menor ángulo posible para obtener la mayor dd 


sensibilidad. 0<9<20" 
0<0<20" "| ma 


El cómputo de la velocidad se realiza con la expresión ya o. 2 80  Blgrados] 


vista 
772 ESC 
2 f cosO 


Si mantenemos constante c, f,, cos € y definimos una constante como cte = c/2 f, cos 6, 


entonces, 
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PW direccional: Esquema en Bloques y Operación (1) 


Lave 


> a IP Demodulador > B 


Básicamente se trata del sistema no direccional al que se ha agregado un demodulador, una 
unidad de muestreo y retención y un filtro pasa banda (antialiasing). Las unidades de 
muestreo y retención reemplazan a la llave L, . Los filtros de hasta 8 polos están 
encuadrados en una ventana frecuencial comprendida entre 100 y 800 Hz. Estos sistemas 
posen gran bW e alto slew rate. 


Las señales provenientes de cuerpos en estado estacionario y de pequeñas señales de los 
mismos en movimiento son amplificadas con la misma ganancia, el receptor requiere de un 
gran rango dinámico para evitar saturación de señales provenientes de las paredes de las 
arterias. Por ello, el receptor debe tener también la capacidad de recuperarse rápidamente de 
la saturación y disponer de un ancho de banda que es inversamente proporcional a la 
resolución axial del sistema. 
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Limitaciones del PW Direccional (1) 


El Doppler pulsado presenta dos grandes limitaciones: en primer lugar, igual que el CW la 
desviación Doppler es proporcional al coseno del ángulo que forma el haz ultrasónico con el 
flujo sanguíneo. Si el haz es paralelo al flujo, el ángulo es cero y el coseno 1, de manera que 
la desviación Doppler representa verdaderamente la velocidad del flujo. Con otros ángulos, 
esta premisa ya no es cierta. En la práctica no se utilizan ángulos mayores de 20* ya que 
producen estimaciones del flujo significativamente más bajos de lo que son en realidad. 


Ejemplo 

Estimar el error cometido en una medición cuando el ángulo que forma el haz ultrasónico y 
el flujo es de 2, 5, 20 30 y 45 grados. El instrumento está calibrado para la condición ideal, 
8 = O grados. Para esta condición, V = 0.5002 m/s, f¿ = 3.126 kHz, f, = 5 MHz y c = 1.600 
m/s. Por lo visto, se tiene: 


V,, =0,5005, V,, =0,5021, V,, =0,5323, V,. =0,5775, V,., =0,707 


En función de estas cifras, se deja que el alumno estime el error cometido para 6 > 0. 

En segundo lugar, la máxima velocidad del flujo que puede medirse viene determinada por la 
frecuencia con la que el sistema emite un pulso y abre la ventana de medición (t, en la Fig. y 
por la profundidad a la que se encuentra el volumen de muestra. A frecuencias bajas, el 
volumen de muestra es pequeño, por lo que, la posibilidad de encontrar anormalidades en la 
velocidad y dirección de los eritrocitos es escasa. A frecuencias mayores puede ocurrir que un 
segundo impulso sea emitido antes de que el eco del primero haya llegado de retorno al 
transductor. Por otra parte, cuando mayor sea la profundidad a la que se encuentra la 
muestra, más tiempo tarda el eco en volver, por lo que puede ocurrir que, la señal de retorno 
sea interferida por un nuevo pulso antes de que pueda ser captada por el transductor. 
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Limitaciones del PW Direccional (2) 


En la práctica, la máxima velocidad del flujo sanguíneo que puede ser determinada sin 
ambigúedad, viene dada por el límite de Nyquist que está dado por la frecuencia de repetición 
de pulso y la profundidad de la muestra. Si este límite se supera, se produce una saturación de 
la señal o "aliasing". El aliasing ocurre a partir de los 0.40-0.60 m/s [4,6,7], siendo así que el 
flujo cardiaco normal puede alcanzar 1.5 m/seg y en algunas patologías hasta 6 m/s. En estos 
casos el PW no tiene utilidad y hay que recurrir al CW 


Ambigúedad en la velocidad y en el rango 
Desplazamiento Doppler mad comparar la relación de fase con la señal de referencia OM. 
El máximo cambio de fase que puede obtenerse entre dos pulsos está limitado al intervalo -TT a 


+TT radianes (las medidas angulares se repiten cada 211 radianes). Por lo tanto la frecuencia de 
muestreo debe ser por lo menos el doble de la componente de frecuencia más alta. 


donde: fama, es el máximo desplazamiento Doppler que puede ser detectado 
Pomo = di sin ambigúedad y f, la frecuencia de repetición del pulso o frecuencia de 
SN muestreo. V..., sin ambiguedad, es: 
ds E fo = F. 2f Va, cosO 
d max A O O 


2 C 


La posición de la llave L, está determinada por el tiempo de demora, entre la transmisión y el 
comienzo del tiempo de adquisición de la señal recibida. Hay ambigúedad en la determinación 
del rango debido a que las señales que arriban al transductor, pueden ser ecos resultantes del 
pulso de transmisión anterior, localizándose alrededor de: cl 2d ] 
dn = == 
2 


donde: n= 0 o entero positivo y ts tiene que ver con PRF 20 


Limitaciones del PW Direccional (3) 


Haciendo los reemplazos correspondientes, Amax Y Vmax a detectar sin ambigúedad: 


se trata de un sistema 


C 
d max == acoplado, donde se debe 
2 > ms 
: Se ti a lograr una situación de 
cf e tiene que — ' Y 
e s 8 f.cose| compromiso entre Vmax 
4 f,cosO y la profundidad máxima 
de trabajo. 
Vmax 
[mm/s] 
6 
En la práctica, los fabricantes de 
equipos los diseñan de acuerdo a 4 Le 
gráficas similares a ma), 
Es decir, Vmax = f(d) con f, como > NUS 
parámetro 5 
, 4 8 12 Distancia 
Ejemplo [mm] 


A partir de los datos obtenidos en el ejemplo anterior, estimar la d,,.. Y V,, para que no 
exista ambigúedad si f, = 5 MHz. Remplazando valores: 


fama = 17,24KHZz 
max =23,2cm 


Vox = 2,8m/s 2 


Resumen de las diferencias entre CW y PW 


l  Doppler Continuo (CW) | Doppler Pulsado (PW) 


Brinda distribución de velocidades de | Brinda distribución de velocidades 
Ventajas los glóbulos rojos de los glóbulos rojos 
Permite determinar velocidades Permite determinar velocidades 
máximas, mínimas promedio máximas, mínimas y promedio 
Presenta variación temporal del 
flujo 
Mejores posibilidades de acceso a El volumen de la muestra puede 
zonas pequeñas definirse a diferentes 
profundidades 
Permite estudiar sitios 


Emite y recibe continuamente 
determinados: resolución de rango 
Puede medir altas velocidades 9 


No permite estudiar sitios determinados: 


ambigúedad de rango ecocardiográfica para ubicar su 


posición 
Emite y espera el regreso del eco 
para hacerlo nuevamente 
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La velocidad máxima está limitada 
Desventajas por la profundidad: se produce 
aliasinc 
Siempre requiere información 


Algunos Ejemplos 


A la izquierda un eco doppler normal, a la derecha, un eco doppler anormal. Se 
observa, por un lado, el ensanchamiento de la base del registro y por otro, el 
incremento del reflujo 


Ejemplos de Eco doppler 


A A LL A A Md A 
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+” 


PLA IE. ARAFAT 


uu 
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Otras Limitaciones del PW Direccional 


La máxima velocidad del flujo que puede medirse viene determinada por la frecuencia 
con la que el sistema emite un pulso y abre la ventana de medición (td) y por la 
profundidad a la que se encuentra el volumen de muestra. 


A frecuencias bajas, el volumen de muestra es pequeño por lo que la posibilidad de 
encontrar anormalidades en la velocidad y dirección de los eritrocitos es escasa. 


A frecuencias mayores puede ocurrir que un segundo impulso sea emitido antes de 
que el eco del primero haya llegado de retorno al transductor. 


Por otra parte, cuando mayor sea la profundidad a la que se encuentra la muestra, 
más tiempo tarda el eco en volver, por lo que puede ocurrir que la señal de retorno sea 
interferida por un nuevo pulso antes de que pueda ser captada por el transductor. 


En la práctica, la máxima velocidad del flujo sanguíneo que puede ser determinada sin 
ambigúedad viene dada por el límite de Nyquist determinado por la frecuencia de 
repetición de pulso y la profundidad de la muestra. Si este límite se supera se produce 
una saturación de la señal o "aliasing". 
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Variantes de PW (1) 


Existen tres variantes básicas: 
* Doppler Color 

* Doppler de Energía 

* Doppler de Pared 


El Doppler Color (DC) utiliza el efecto Doppler 
generado por los glóbulos rojos en movimiento, 
utiliza el color para codificar la frecuencia de 
corrimiento Doppler. El color y los diferentes tonos 
de los mismos son determinados por la dirección y 
la magnitud del cambio de frecuencia. AMPLITUD 
Para codificar la turbulencia se usa la varianza, la 
que expresa las desviaciones de las diferentes 
velocidades presentes con respecto a la media. Y VARIANZA 
Como generalmente la profundidad y la resolución 
de velocidad son los parámetros más importantes, 
éstos se mantienen sacrificando el ángulo X 
sectorial, la frecuencia de imagen y la densidad de 

líneas. Los mapas en color del flujo se superponen 

sobre los ecocardiogramas que se muestran en 

gris. 


POTENCIA 
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Variantes de PW (2) 


El Doppler de Potencia no tiene problemas de aliasing ni dependencia angular, 
permite determinar el flujo de pequeños vasos y de baja velocidad. Recibe 
diferentes nombres: Color Power Doppler, Ultrasound Angio, Color Doppler Energy, 
Color Power Angio. Asigna valores a diferentes amplitudes de la señal Doppler. La 
amplitud está determinada por la concentración de los dispersores (eritrocitos) que 
producen el desplazamiento de frecuencia. La amplitud está relacionada con el 
volumen en movimiento. Integra el área bajo el espectro. La ventaja respecto a 
otros modos es la presentación uniforme del flujo, aún en aquellas zonas de 
incidencia normal. Sin embargo se pierde información de dirección, velocidad y 
característica del flujo. 


Distribución de Potencia Pseudo-Instantánea 


Potencia 


D.3 
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Frecuencia (Hz) Tiempo (s) 


Variantes de PW (3) 


Brinda distribución 
Ventajas de las velocidades espacial de las partículas 
en movimiento 
Permite calcular para cada 
ubicación el valor medio, la flujo 
varianza y el signo de la 
velocidad ángulo 


Desventajas No realiza una No permite evaluar 
caracterización completa | ninguna característica del 
del flujo flujo 


El Doppler Pared realiza el análisis del movimiento de las paredes cardíacas, 
donde se tiene baja velocidad y alta amplitud del flujo, a diferencia del resto del 
organismo donde se tiene alta velocidad y baja amplitud. Debido a que se trata de 
flujo turbulento y rotacional, el sistema dispone de filtros de amplitud para suprimir 
las señales producidas por el flujo, en general se tratan de filtros adaptivos. 
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Ejemplos de Combinación Doppler-Imagen 


PRF = 1.5 KHz 


- 19.2 
ecmís 
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Introducción 


A diferencia de las áreas de los END y del diagnostico médico, los transductores para 
potencia necesitan de un tratamiento diferente. Los primeros, están constituidos por 
una cerámica, un material reflector (baking), un circuito para adaptación de 
impedancias y en algunos casos de una línea de demora o zapata para realizar 
incidencia angular. En el caso de los transductores para potencia, la realidad es un 
poco diferente. Por un lado, se necesitan al menos dos cerámicas y a continuación de 
éstas se deben acoplar dos masas metálicas, que en general son de materiales y 
dimensiones diferentes. Para que el transductor sea operativo en lo que respecta a la 
máxima potencia entregada al medio, se necesita por un lado, disponer de un circuito 
pasivo para la adaptación de impedancia entre éste y el generador de señal y por 
otro, de un tornillo o perno que mantiene al conjunto cerámica-masas unido con una 
tensión mecánica específica. 

Las cerámicas piezoeléctricas presentan mejores factores de conversión 
electromecánica que los materiales para transductores magnetoestrictivos. En estos 
últimos, los campos alternos dan lugar a corrientes parásitas, que producen 

pérdidas muy dependientes de la frecuencia y de las corrientes parásitas inducidas 
en el material y que son de considerable magnitud. Por ello, los transductores 
magnetostrictivos son dispositivos especialmente orientados para trabajar en baja 
frecuencia entre 20 y 25 kHz. Los componentes metálicos usados en la construcción 
del transductor modifican y afectan a la capacidad de potencia y rendimiento del 
mismo. En general el empleo del aluminio, duraluminio y acero son los materiales 
metálicos que más se usan en los transductores de potencia. 232 


Introducción 


Las aplicaciones del ultrasonido se clasifican generalmente en dos grandes grupos 
que se refieren, respectivamente, al uso de bajas y altas intensidades. La mayoría de 
las aplicaciones del ultrasonido de potencia se llevan a cabo en el campo de 
frecuencia entre 18 y 200 kHz y con intensidades acústicas que varían desde 
aproximadamente 0,1 W/cm* hasta varios KW/cm? permitiendo obtener altas 
aceleraciones y concentración de energía en medios de muy diversa naturaleza. 

Se los puede clasificar en dos áreas perfectamente definidas: 


Los que producen modificaciones permanentes en el medio sobre el que se 
actúa. Existe una amplia gama de procesos que pueden ser provocados mediante 
ultrasonido de potencia. Los más importantes en el área industrial son: mecanización, 
soldadura, formación de metales, etc. en medios sólidos; limpieza, aceleración de 
reacciones químicas, emulsificación, atomización, desespumación, secado, 
aglomeración de aerosoles, etc. en medios fluidos. En el área de la medicina se 
emplean en la eliminación de cálculos como es el caso de la litotricia. 


Los que no producen modificación alguna en el medio sobre el que se actúa. El 
caso típico de esta aplicación es la del SONAR. La energía acústica enviada al medio 
se emplea para determinar si en su trayecto existe algún elemento que produzca un 
eco. Se emplea en la detección de banco de peces y como ayuda a la navegación. La 
potencia del sistema se utiliza para compensar la atenuación debida a la distancia y a 


la reflexión sobre el blanco. 
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Materiales para Transductores de Potencia 


1. Cerámicas piezoeléctricas 
2. Materiales piezomagnéticos 
Diferencia ente materiales piezoeléctricos y piezomagnéticos 


Cerámica Niquel- Terfenol 
piezoeléctrica cobalto 
Factor de acoplamiento 0,7 0,5 0,7 
Máxima deformación 200 39 1000 
(ppm) 
Módulo de Young (Gpa) 74 190 26 
Densidad (g/m? 7.0 8,9 9 
Componentes metálicos 
Material Módulo de Densidad Velocidad Resistencia Máxima 
Young 10%Kg/m? del Sonido ala Fatiga deformación 
(1011N/m2) (m/s) (107N/m2) (103) 
Aleación de 
Ti 1,06 4,4 4900 72 6,8 
(Ti-6A1-4V) 
Duraluminio 0,74 2,8 5130 19 2,6 
Aluminio- 1,43 8,5 4070 37 2,6 
bronce 
Latón 0,89 8,4 3240 15 1,5 
Acero de 2,2 7,8 5200 55 210 
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Variantes de Transductores de Potencia (1) 


Ejemplo de Transductores Electromagnéticos 


Radiación 
Eléctrico Mecánico acústica 


S 
RR o” E: * 
p> Br 
nl 2 
j | 
e¿Bl 1 
| (RR FL 
0 e 


se... 


Y, 


Y 


Br 


(R¿+Rg) 


2. «Suspensión de borde flexible 
' 3, Suspensión central flexible 


ña .. Bl 
| 1, Sonodeflector infinito sap! 


4, Bobina móvil 
S. Conexiones eléctricas (b) 
' 6, Orificios de salida de aire, 
| 7, Cono (diafragma! 
3. Imán permanente 
1 9 Estructura de soporte Espesor del borde 0,25 mm 


Borde de fijación 
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Variantes de Transductores de Potencia (2) 


Transductor Piezoeléctrico para Sonar 


compuesto 


perspex de moldeo 


conductor 
eléctrico 


fondo de 
goma 


conductor 
eléctrico 


coaxil 
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Transductor Langevin (1) 


El transductor piezoeléctrico más característico empleado en aplicaciones de potencia 
es el que se conoce como transductor Langevin, cuya estructura proviene de la de un 
transductor piezoeléctrico de cuarzo diseñado originalmente por Langevin para 
producir cavitación. Cuando se emplearon por primera vez las cerámicas 
piezoeléctricas para traductores ultrasónicos de banda estrecha a bajas frecuencias, el 
transductor consistía en un simple bloque cerámico. Sin embargo, esta solución no 
resultó ser muy útil debido a la baja resistencia mecánica, a la tensión de las 
cerámicas y a las dimensiones necesarias para que el bloque cerámico fuese 
resonante a bajas frecuencias. Este transductor, consta de una estructura formada por 
dos anillos cerámicos del mismo espesor, dispuestos entre dos cilindros metálicos, 
formando una estructura que vibra extensionalmente a media longitud de onda. Esta 
estructura queda fijada a través de un tornillo que aumenta la resistencia a la tensión 
de las cerámicas y permite trabajar a los distintos materiales como si fueran uno. 


contramasa a cabeza 
«——— tornillo 
| Ejemplo de transductores Langevin 
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Transductor Langevin (2) 


Comportamiento de un par de discos piezoeléctricos conectados en serie con la 
polaridad longitudinal invertida. Las flechas representan la dirección, de la polarización 
de los discos piezoeléctricos 
Partes Constitutivas 


Contramasa (Backing): Cilindro metálico 
posicionado en la parte trasera del 
transductor. El material debe poser alta 
densidad y actúa como un gran atenuador, V = Vmaz V= V = —Vinar 
absorbiendo /a energía que se produce en 
la parte trasera del mismo. De la 
adaptación de impedancias entre la 
contramasa y el elemento activo 
dependerá la resolución y la sensibilidad 
del transductor. 


Cabeza (Matching): Cilindro metálico 
ubicado entre el elemento activo y el 
medio de transmisión, cuya función 
es de servir de transformador 
acústico acoplando las impedancias 
acústicas del elemento activo y el 


Buscando mejorar el acoplamiento medio. El espesor necesario para 
acústico entre las cerámicas y los conseguir este efecto es de 1/4 de la 
bloques metalicos y prevenir fallas por longitud de onda emitida. El material 
fatiga mecánica, se ha visto la necesidad más común para este tipo de placa es 
de aplicar una presión entre 89 y 97 MPa el aluminio con una Z de 17,1*108 
para mantener unidas las masas metálicas Kg/m?s. 


con los cristales. Para ello se usa un 
tornillo central. 238 


Transductor Langevin (3) 


Análisis detallado de un transductor tipo Langevin 


Como se desprende de lo visto, la base del transductor Langevin tiene como punto de 
partida el análisis de una barra sometida a esfuerzos cuando se la hace vibrar . donde 
se ponen de manifiesto las impedancias puestas en juego. 


Impedancias 
Za = Zaire 
Zc =Z cerámica 


Za Zem Zca Zm 


Zce = Z cabeza 

Zem = Z contranmasa 
Se tiene: Zm = Z medio 
Zcm >>> Za, Zem > Zc, Zca >= Ze 


entonces: 
Interfaz Zc — Zem y Zc - Zca ma) Interfaces Duras 
Interfaz Zcm - Za y Zca - Zm mm) Interfaces Blandas 


En la interfaz Zcm - Za existirá reflexión total y por lo tanto inversión de fase de 
la onda reflejada 
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Transductor Langevin (4) 


2/4 


Acción de la contramasa 


Las componentes reactivas han de cancelarse 
mutuamente. Entonces, el sistema vibrante 
contramasa-cerámica se comporta como una lámina 
A/4. Cuando se excita al transductor, parte de la 
vibración producida se dirige hacia la cabeza radiante 
y parte a la contramasa. Como la longitud de la 
contramasa es de 2/4, la señal llega a la interfaz 
contramasa-aire con una diferencia de fase de rr/2. En 
ésta se produce reflexión total. La onda reflejada 
llegará a la interfaz cerámica-cabeza radiante con una 
diferencia de fase de 3/2rr. La radiación hacia 
adelante se compone en fase con la reflejada en la 
cara trasera del transductor, incrementando de este 
modo, la radiación al medio. 


“Sándwich” 


2/2 


1 


31/4 


3IT 


a 
na 
1. 
A 


Xx Us 


IO A AS MM. 


72 


sd 
Aa 
E 
mn e 
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Transductor Langevin (5) 


Relacionando el transductor con ondas longitudinales en una barra, las impedancias 
contramasa-cerámica y cerámica-contramasa están dadas por: 
k es el número de onda, | espesor del medio 
== c. te “(kl >) considerado. Ambas impedancias son reactivas 
cm-ce IP cm 8 cm 
puras y se cancelan para un espesor de 2/4 


a Si =) ¡E Cee 18 (kl) 


Dividiendo Zmc-ce por Zce-cm y multiplicando numerador y denominador por el área 
de cada elemento: 


y 0 Pon Com 


= te(kl )te(kl >) (1) Ecuación de Langevin 
o Pee Cos Ñ ll 


La determinación de la longitud de la cabeza y contramasa del transductor está 
basada en el conocimiento de la longitud total de la oscilación, que debe ser de M2. 
En el caso de metales se la puede expresar como 


donde: c, es la velocidad del sonido en una barra de metal 
de diámetro d y longitud |, fr la frecuencia de resonancia 
del transductor y A = v,/fr. 241 


Transductor Langevin (6) 


La longitud de un transductor simétrico se puede expresar como: 


l=21 +21. (2) 
|. es la longitud de una de las masas, por ejemplo, la contramasa, entonces: 


l 
n= 0-21.) 


A es una constante que depende del tipo de metal y de las cerámicas usadas, como 
también de la frecuencia de trabajo. Si en la expresión anterior se sustituye l' de la 
expresión (2): 

—= C, 


“Ah 
Válida para | > d. En la práctica el valor de la cte. A se puede tornar un poco 


complicado, por ello, y como primera aproximación suele darle el valor de 0,5, 
quedando la Ecu. (2): 


| = e 
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Transductor Langevin (7) 


Ejemplo 

Se dispone de dos cerámicas con las siguientes características: p = 7500 kg/mY, 
módulo de Young (E) = 7,4*101% N/m?, diámetro externo 38,1 mm, diámetro interno de 
12,7 mm espesor 6.35 mm, frecuencia de trabajo f = 22,5 kHz. El material de la 
cabeza y contramasa es acero inoxidable con un diámetro de 40 mm y p = 7715 
kg/m*, módulo de Young (E) = 2,18*10*! N/m?. Determinar las dimensiones del 
transductor simétrico, primero de acuerdo a la expresión (1) y luego con la (2). 


*3(y.0 31011 
ces O 10 = 3141 m/s AN as 10 = 5315 m/s 
P 7500 T715 


(40-10) 7/4 


z Zen" Pen" Con" Mem 1715-5315 Ar _ 17406 E dai 
ZA. p.oc-:A. 7500-3141 A, (38,1-10? -12,7-107*) 7/4 
Reemplazando Z en la (2) lo. = Ly = arctan g| —= po 
j ¡ "Lag. LI] Kon 
Ica = lcm= arctan g sa 5 E j - =56,5-107m=56,5mm 
te (27 -22,5-10*-6,35-107/3141-) | 27 -22,5:10*/5315 


lla = lo = 56,5 mm 


Longitud total del transductor = 2*56.5 + 2*6,35 = 125,7 mm 243 


Transductor Langevin (8) 


A partir de la ecuación (3) 


| ==" -051 
4f 
C.,= 5315 m/s 
f, = 22,5 KHz 
L., = 6,35 mm 
5315 


O 


Longitud total del transductor = 2*55.9 + 2*6,35 = 124,5 mm 


Al principio del proyecto, el error cometido con esta aproximación es aceptable. 
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Transductor Langevin (9) 
Transductor + Amplificador Mecánico 


Para este tipo de amplificador, la relación entre las impedancias y las velocidades está 
dada por: 
2/4 1/4 


Y _Z,_58,Y,c fg(w,/2c,) 
A VW» Z, SY c,tg(w,/2c,) 


AL2 


El factor de amplificación es la relación de áreas de los extremos del amplificador. 
Si la relación D,/D, = 2, la relación v,/v, = 4 


Transductor Langevin multicapa y terminado con amplificador mecánico 


CM CA AM 


A/2 


Ejemplos de modelos de 
amplificadores mecánicos » 


3 
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Transductor Langevin (10) 


Ejemplos de simulación por elementos finitos de la deformación 
mecánica del transductor excitado 
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Transductor Langevin (11) 
Introducción al modelado eléctrico de un transductor Langevin 
Modelo de Mason de un transductor piezoeléctrico 
Co Cu -Cm Ce L 


%1 


HH 
. L 


E y C, son los elementos eléctricos del sistema, C;, L y R los mecánicos y 
Cy el electromecánico. E es la diferencia de potencial aplicada a los electrodos 
de entrada del transductor. C, representa la capacitancia electrostática del 
capacitor formado por los electrodos y el material piezoeléctrico, todo asociado a 
la permitividad del material e,, que en el caso unidimensional se conoce como 
constante dieléctrica y e, la permitividad del vacío 


C Ey E, A f= 4p,!, A  Lesel equivalente mecánico de la masa de la 
07 — 2 cerámica, asociado a la energía inercial mecánica 
le T 
Os S l, Cp es la capacidad de almacenamiento de energía potencial 
E A mecánica del material 247 


Transductor Langevin (12) 


Introducción al modelado eléctrico de un transductor Langevin (2) 
Modelo de Mason de un transductor piezoeléctrico 


C. — Er E, $ Es la razón de carga eléctrica aplicada a la cerámica y la fuerza 
M= d requerida para que éste no modifique sus dimensiones, 
contrarrestando el efecto piezoeléctrico 


Constantes Unidades: 
| desplazamiento dieléctrico desarrollado | coulomb / metro? C 
tensión mecánicaaplocada Pa > N 
E = cons tan te 
tensión mecánica desarrollada metros / metro 
campo eléctricoaplocado volts / metro 


T = cons tan te 
a 


campo eléctrico desarrollado Volts | metro 


tensión mecánica aplocada Pa 
D = cons tan te 
tensón desarrollada metros / metro 
desplazamiento aplicado al dieléctrico | coulombs / metro” 


T' =cons tan te 
Tabla:3.1. -Definición-de-las-constantes-d-y-gf 
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Transductor Langevin (13) 


Introducción al modelado eléctrico de un transductor Langevin (2) 

Modelo de Mason de un transductor piezoeléctrico 

Ejemplo 

Mediante la simulación del modelo Mason, obtener la respuesta en frecuencia de un 
transductor tipo PZT que tiene los siguientes parámetros: ey = 1750, e, = 8,85*10112, 
p = 7500 Kg/m?, O = 0,3, kp = 0,55, kzz = 0.64, s = 7,4*10*0, diámetro interno (di) = 
12,7 mm, diámetro externo (de) = 38,1 mm, espesor (lp) = 6.35 mm, d = 1,2*10"*? y 
R =1,7 KO. 


Cálculo del área Co = 9,88 nF 
L=78,2mH E 9 


Modelo de simulación 


Probe 1-NODE 


1.7k 
R1 
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El principio del sonar se basa en la técnica de pulso-eco. Por razones de economía 
se usa un solo transductor para transmitir y recibir, tal como se vio a lo largo de esta 
obra. A diferencia de la técnica pulso-eco conocida para aplicaciones en END, al 
transductor no se lo excita con un pulso angosto sino con un ciclo (o varios) de una 
senoide, similar a los sistemas de Doppler pulsado. El generador de pulsos controla 
el periodo entre los mismos y su ancho que dependerá de la aplicación específica. El 
receptor que es similar a los utilizados en pulso-eco, está formado por un 
amplificador pasa banda sintonizado a la frecuencia de trabajo del transductor, un 
detector de envolvente y un amplificador final necesario para llevar la amplitud de la 
señal detectada a los niveles necesarios para el procesamiento e indicación. 


Esquema en bloques de un sistema para sonar 


Oscilador 
Amplificador 
Sintonizado 


Generador de 
Pulsos 


Amplificador 


Sistema 
Indicador 
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Sonar (2) 


Las altas frecuencias son muy adecuadas cuando se desea enviar pulsos cortos, 
que son de particular interés en aguas poco profundas, debido a que tanto la 
distancia mínima de medición y resolución dependen de la duración del pulso de 
ultrasonido. Por otro lado, la absorción del sonido en el agua de mar aumenta 
abruptamente con la frecuencia, por lo tanto debe encontrarse el mejor compromiso 
entre ambas. Para un sonar con rango máximo de 100 m, la frecuencia óptima se 
encuentra entre 150 y 200 kHz. La distancia mínima que se puede medir es 
típicamente 30 cm; con una resolución de aproximadamente la mitad de esta cifra. 
A medida que los pulsos de ultrasonido se hacen más cortos, el bW del transductor 
tiende a transformarse en un transductor de banda ancha. Si esto sucede, la 
selectividad y la SNR se verán afectadas. El bW requerido a -3 dB es de 
aproximadamente 10 a 15 kHz correspondiente a un factor de calidad mecánica de 
Q = 15. 

Los transductores de flexión o a diafragma (electrostáticos) se prefieren para 
aplicaciones de largo alcance operando a frecuencias inferiores a 50 kHz. La 
frecuencia de resonancia se rige principalmente por el espesor del disco, pero la 
directividad se rige por su diámetro, la variación de la amplitud de la vibración a 
través de su superficie y la frecuencia de trabajo. 
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Sonar (3) 


Los transductores piezoeléctricos como los vistos a lo largo del curso con algunas 
modificaciones se usan para sonar 


compuesto 
de moldeo 


perspex 
Esquema de un transductor RC 
piezoeléctrico basado en PZT eléctrico 


fondo de 
goma 


conductor 
eléctrico 


cable 
coaxil 


Al seleccionar el modelo de transductor, las siguientes propiedades deben ser 
consideradas: 

- Baja absorción de sonido. 

- Poder adhesivo y la elasticidad (protección contra daños mecánicos). 

- Tolerancia al agua de mar y luz solar. 


- Superficie lisa para reducir el crecimiento de las algas y facilitar la limpieza. 259 


Parámetros de un transductor para sonar típico 


Frecuencia de operación f, 151 kHZ 


Capacitancia del disco a 1 kHz 640 pF 


Inductancia de compensación L,, Depende de la capacitancia del disco más la del cable 


Impedancia a f¿ (con Lp) 1,13K0 


bW a -6 dB (sin resistor de carga) 192 (diagrama de radiación) 


bW a -6 dB (con resistor de carga) 132 (diagrama de radiación) 


Sensibilidad (eco a 27 cm) 0,07 
tores para sonar Murata 


PW mínima del pulso de excitación 80us * 12 cm 


Ejemplo de transductores comerciales para sonar (Firma Murata) 


Frecuencia Capacitancia Impedancia Directividad Potencia de 
(resonancia) (13) (ES TELE) (grados) entrada 
(kHz (69) (w) 


UT200BA8 200 1700 310-590 22 50 


UT200LF8 200 2700 230 - 430 12 200 
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